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CONVERSION  OF  MICOM,  TIME-PHASED,  LIFE-CYCLE, 

COST-ESTIMATING  MODEL  FROM  COBOL  TO  FORTRAN  IV 

ABSTRACT 

LICEM  is  a  computer  model  which  may  be  used  to  generate  Time- 
Phased  Life-Cycle  Cost  Estimates  (LCCE)  for  personnel  or  materiel 
systems.  The  input  to  this  model  is  in  a  form  compatible  with  the  Army 
Materiel  Command's  Improved  Cost  Estimating  Project,  Phase  III  (ICE-III). 

The  cost  for  a  system  can  be  computed  for  as  many  as  thirty 
equal  increments  of  time,  and  can  be  summarized  in  up  to  nine  levels  of 
complexity.  The  model  estimates  a  system  cost  for  each  time  increment 
as  well  as  the  total  cost  over  the  life  of  the  system. 

The  model  is  written  in  FORTRAN  IV  specifically  for  the 
Ballistic  Research  Laboratories'  Electronic  Scientific  Computer 
(BRLESC) . 
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CONVERSION  OF  MICOM,  TIME-PHASED,  LIFE-CYCLE, 
COST-ESTIMATING  MODEL  FROM  COBOL  TO  FORTRAN  IV 


1 .  INTRODUCTION 

The  Cost  and  Analysis  Division  of  the  U.S.  Army  Missile  Command 
(MICOM)  developed  a  Time-Phased,  Life-Cycle,  Cost-Estimating  Model  which 
has  gained  wide  acceptance  throughout  the  Army  Materiel  Command  (AMC) 
and  its  subordinate  commands  for  assessing  the  life  cycle  costs  of  Army 
Systems  (Reference  1).  The  model  is  especially  useful  in  that  the  input 
format  corresponds  with  that  described  in  AMC's  Project  ICE  (Improved 
Cost  Estimating)  Phase  III. 

The  program  was  originally  written  in  the  COBOL  computer 
language.  Because  of  a  need  for  this  model  in  the  Systems  Methodology 
and  Resource  Studies  Office  (SM^RSO)  of  the  U.S.  Army  Materiel  Systems 
Analysis  Agency  (USAMSAA)  and  the  lack  of  COBOL  capability  on  the 
Aberdeen  Research  and  Development  Center's  (ARDC)  computer,  the  model 
was  converted  to  FORTRAN  IV. 

This  memorandum  describes  the  FORTRAN  IV  version  of  MICOM 's  Time- 
Phased,  Life-Cycle,  Cost-Estimating  Model. 

There  are  many  cost  measures  which  can  be  associated  with  an 
Army  weapon/support  system.  Examples  are  procurement  costs,  operating 
costs,  life-cycle  costs,  etc.  The  life-cycle  cost  estimate,  associated 
with  each  system,  describes  the  cost  of  that  system  from  its  R$D 
status  through  its  operational  phase  and  to  its  retirement,  and  there¬ 
fore  provides  the  most  comprehensive  portrayal  of  the  cost  of  the 
system. 


Realizing  that  life-cycle  costs  were  playing  an  important  role 
in  its  acquisition  process,  the  Army  recognized  the  need  to  standardize 
the  procedure  for  determining  these  life-cycle  costs.  This  procedure 
was  defined  in  the  Improved  Cost  Estimating  Study  (ICE-III) (Reference  2). 


2.  IMPROVED  COST  ESTIMATING  (ICE-III) 

Army  Regulation  37-18  "...establishes  a  set  of  cost  categories 
and  elements  to  be  used  by  weapon/support  system  cost  analysis 
activities...  These  categories  and  elements  will  be  used  for  life-cycle 


^U.S.  Army  Missile  Command,  Time  Phased  Life  Cycle  Cost  Model,  Cost 
Analysis  Division,  30  June  1969. 

2 

Letter,  Office  of  AMC  Comptroller,  Program  for  Improved  Cost  Estimating 
(ICE)  Phase  3,  14  May  1970. 
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analyses  of  existing  and  proposed  weapon/support  systems,  to  establish  a 
uniform  basis  for  cost  analysis...  and  are  to  be  used  in  life  cycle 
costs  analyses  for  cost  effectiveness  and  decision  oriented  studies... 

The  objective  of  this  regulation  is  to  improve  cost...  estimating  proce¬ 
dures  within  the  Department  of  Army  through  the  use  of  uniform  categories 
and  elements  for  weapon/support  system  costs."  (Reference  3). 


AR  37-18  defines  Cost  Categories  as:  "The  major  divisions  of 


weapon/support  systems'  cost  from  inception  to  retirement  of  the 
system."  (Reference  3).  The  following  cost  categories  are  defined  in 
detail  in  AR  37-18: 


•  Research  and  Development 

•  Investment  Non-Recurring 

•  Investment  Recurring 

•  Operating 

Cost  Elements  are  defined  as  "The  subdivision  of  cost 


categories  related  to  work  areas  or  processes  performed  in  developing, 
producing,  and  operating  a  weapon/support  system."  (Reference  3). 


A  complete  breakdown  of  the  cost  structure  associated  with 


weapon/support  systems  is  given  later  in  this  report .  The  cost  struc¬ 
ture  proposed  in  AR  37-18  forms  the  basis  for  the  rationale  used  in  this 
cost  estimating  model 


3.  DESCRIPTION  OF  ICE-II I  CODES  AND  INDEXING  STRUCTURE 


The  chart  on  the  following  page  describes  the  ICE-II I  codes  and 


indexing  structure.  These  are  used  to  create  codes  for  each  level  of 
data  that  are  input  to  the  life-cycle  cost  model.  A  summary  of  these 
codes  is  now  given: 


Pair  of  Digits 


Represents 


1 


The  system  under  consideration 


4 


2 


3 


Cost  Category 
Cost  Element 
Type  of  Cost 


7 

Army  Regulation  AR  37-18,  "Weapon/Support  Systems  Cost  Categories  and 
Elements,"  2  July  1968. 
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ICE-PHASE  III  AND  AR  37-18  CODES  AND  INDEXING  STRUCTURE 
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7-9 


Sub -Element 
Appropriation 

Work  Breakdown  Structure  of  System 


The  code  associated  with  each  level  of  data  contains  N  pairs 
of  digits  where  N  9.  For  illustrative  purposes  we  will  now  describe 
levels  of  data,  and  assign  their  associated  level  codes  from  the  index¬ 
ing  structure. 

Example  1:  Suppose  a  level  in  which  we  are  interested  is  the 
RDT§E  costs  of  direct  labor,  of  in-house  engineering  in  the  development 
cost  category.  The  associated  code  would  be  as  follows: 


Code:  01 

01  01  02  01  01 

Pair 

No. :  1 

2  3  4  5  6 

Pair  1: 

System  code 

Pair  2: 

Development  cost  category 

Pair  3: 

Engineering-cost  element 

Pair  4: 

Type  of  cost-in-house 

Pair  5: 

Sub-element-direct  labor 

Pair  6: 

Appropriation-RDT§E  costs 

Example  2:  The  level 
under  contract  for  crew  initial 

code  for  the  MPA  costs  of  Direct  Labor 
training  is: 

Code:  01 

03  09  01  01  04 

Pair 

No:  123456 


Pair  1: 

System  code 

Pair  2: 

Cost  category-investment  Recurring 

Pair  3: 

Cost  element-initial  training  crew 

Pair  4: 

Type  of  cost-contract 

Pair  5: 

Sub-element-direct  labor 

Pair  6: 

Appropriation-MPA 

In  this  manner  the  codes  are  assigned  to  each  individual  level  of  data. 
Notice  that  no  codes  were  given  for  the  work  breakdown  structure,  since 
it  will  be  peculiar  to  each  system  under  consideration. 
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This  code  and  indexing  structure  allows  a  great  deal  of 
latitude  in  the  amount  of  detail  which  is  to  be  used  in  the  model.  This 
level  of  detail  is  to  be  defined  by  the  individual  analyst  and  will  be 
dependent  upon  his  time  and  resource  constraints  on  data  acquisition. 


4.  THE  MODEL 

The  manner  in  which  the  life-cycle  costs  are  derived  in  LICEM 
is  easily  understood  once  the  cost  breakdown  is  defined  by  ICE-III  and 
AR  37-18.  The  data  are  stratified  in  as  many  as  nine  levels  of  complexity, 
with  the  greatest  complexity  occurring  at  the  work  breakdown  structure 
level,  and  the  least  complexity  occurring  at  the  total  systems  cost  level. 
Since  a  code  is  associated  with  each  level  of  data,  the  final  cost 
associated  with  a  particular  level  is  the  sum  of  the  costs  associated 
with  its  sublevels.  In  this  manner  the  data  are  "summed  up"  until  the 
total  system  cost  is  determined. 

The  costs  for  each  level  of  data  are  entered  into  the  model 
through  a  set  of  input  cards.  There  are  seven  distinct  types  of  cards 
which  could  be  used  to  calculate  the  cost  associated  with  a  level  of 
data.  The  description  of  these  cards,  their  use  and  format,  will  be 
presented  in  Section  9.1,  Preparation  of  Input. 


5.  SUMMING  UP  PROCEDURE 

As  stated  previously,  due  to  the  "summing  up  procedure,"  only 
nine  levels  of  data  are  needed  for  storage  in  core  during  the  operation 
of  the  processing  model.  See  Figure  1  for  visual  description  of  storage 
matrix  used  in  Main  Processing  Program.  The  nine  levels  in  storage,  one 
for  each  of  the  nine  possible  pairs  of  digits  in  the  level  code,  repre¬ 
sent  the  latest  level  of  data  entered,  composed  of  1,2, 3,.. .,9  pairs  of 
digits,  respectively.  As  a  new  level  of  data  (whose  code  is  composed  of 
N(N<9)  pairs  of  digits)  is  entered  and  its  attendant  costs  calculated, 

they  will  replace  the  current  data  in  the  row  of  the  storage  matrix. 
But  before  this  transfer  can  take  place,  the  "summing  up"  procedure 
must  be  performed.  This  procedure  can  best  be  explained  by  example, 
and  one  is  given  in  this  section.  Before  the  example  is  given,  however, 
a  few  words  on  how  the  "summing  up"  procedure  works  in  general  would  be 
useful.  One  can  think  of  all  the  sublevels  of  one  level  (super  level) 
as  being  nested,  or  contained  within  the  super  level.  The  summing  up 
procedure  then  totals  the  contributions  of  all  nested  sublevels,  in  turn 
giving  the  total  for  each  sublevel  in  the  set,  as  well  as  the  total  for 
the  super  level. 
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Figure  1  Pictorial  Description  of  Storage  Matrix 


These  nine  possible  degrees*  of  data  are  all  that  need  to  be 
kept  in  memory  at  any  one  time.  Once  a  level  and  its  sublevels,  if 
any,  have  been  calculated  and  placed  in  the  matrix,  they  are  added  upward 
and  then  written  on  an  output  tape,  thus  clearing  those  rows  for  the 
introduction  of  new  data. 

For  illustration  purposes,  suppose  our  life-cycle  cost  is 
composed  of  five  levels  of  data,  covering  a  period  of  2  years,  that  are 
in  the  first  column  of  Figure  2.  The  times,  noted  along  the  left  margin 
correspond  to  the  entrance  of  the  next  level  of  data.  The  STORAGE  MATRIX 
and  OUTPUT  TAPE  columns  are  snap-shot  views  of  how  the'  matrix  and  tape 
appear  during  each  stage  of  this  processing  example. 

TIME=0:  this  level  is  used  merely  to  show  the  levels  of  data 
that  are  to  be  entered,  and  that  the  istorage  matrix  and  tape  are 
initially  blank. 

TIME=1:  enter  the  first  level  of  data.  Note  that  these  data 
go  directly  to  the  storage  matrix.  Also,  nothing  will  be  placed  on  the 


*The  degree  of  a  level  of  data  is  the  number  of  pairs  of  digits  present 
in  its  level  number,  i.e.,  0102030501  is  of  degree  5. 
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Level  Data  To 

Be  Entered 

Storage 

Matrix 

Output  Tap' 

e 

Code  Costs 

Code  Costs 

01  10  15 

LEVI 

0101  8  9 

LEV  2 

"Blank" 

"Blank" 

TIME=0 

010101  5  6 

LEV3 

010102  4  8 

0102  10  8 

Enter  1st 

LEVI 

01 

10  15 

TIME=1 

Level  of  Data 

LEV  2 

"Blank" 

"Blank" 

LEV3 

"Blank" 

Enter  2nd 

LEVI 

01 

10  15 

TIME=2 

Level  of  Data 

LEV  2 

0101 

8  9 

"Blank" 

LEV3 

"Blank" 

Enter  3rd 

LEVI 

01 

10  15 

TIME=3 

Level  of  Data 

LEV  2 

0101 

8  9 

"Blank" 

LEV3 

010101 

5  6 

Enter  4th 

LEVI 

01 

10  15 

010101  5 

6 

TIME=4 

Level  of  Data 

LEV  2 

0101 

13  15 

LEV3 

010102 

4  8 

Enter  5th 

LEVI 

01 

27  38 

010101  5 

6 

TIMERS 

Level  of  Data 

LEV  2 

0102 

10  8 

010102  4 

8 

LEV3 

"Blank" 

0101  17 

23 

At  this  time  all  data  have  been  entered. 


The  final  summation  procedure  yields 


010101 

5 

6 

Final  Form  of 

010102 

4 

8 

Output  Tape 

0101 

17 

23 

0102 

10 

8 

01 

37 

46 

Figure  2  Level  Data,  Storage  Matrix,  and  Output  Tape  For 
Example  of  Summing  Up  Procedure 
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output  tape  until  one  level  replaces  another  in  the  storage  matrix,  and 
the  "summing  up"  procedure  is  performed. 

TIME=2:  enter  the  second  level  of  data.  This  level  code  0101, 
composed  of  two  pairs  of  digits  is  placed  directly  into  the  second  row 
of  the  storage  matrix. 

TIME=3:  enter  the  third  level  of  data.  This  is  placed  in  the 
third  row  of  the  storage  matrix.  At  this  point  no  "summing  up"  has  been 
done . 


TIME=4:  enter  the  fourth  level  of  data,  010102;  note  that  this 
replaces  010101  in  the  storage  matrix;  hence  010101  must  be  added  up  and 
written  on  the  output  tape.  The  adding  up  is  evident  from  the  values  in 
0101  in  the  storage  matrix,  which  are  now  13,  15.  These  values  are  the 
sum  of  0101' s  original  values  8  and  9  and  the  values  5  and  6  from  010101. 
Finally  010101  and  its  original  values  (5,6)  are  written  on  tape. 

TIME=5:  enter  the  fifth  level  of  data,  0102,  which  will  replace 
0101  in  the  storage  matrix.  Several  things  need  be  done  here. 

•  All  sublevels  of  0101,  (which  is  010102)  must  be  "summed  up" 
and  included  in  0101.  Note  in  the  output  tape  that  this 
yields  a  final  total  for  0101  of  17  and  23. 

•  After  being  summed  up,  these  sublevels  of  0101  must  also  be 
written  on  the  output  tape.  Note  that  level  010102  with 
costs  4  and  8  appear  on  the  image  of  the  output  tape. 

•  Before  0101  is  written  on  tape,  its  final  values  must  be 
summed  to  level  01,  given  01  values  of  27  and  38  in  the 
storage  matrix. 

•  Finally  0102  replaces  0101  in  the  second  level  of  the 
storage  matrix. 

TIME-FINAL:  At  this  time  all  levels  have  been  entered. 

"Summing  up"  the  costs  of  0102  to  01,  giving  totals  of  37  and  46,  and 
writing  these  final  two  levels  onto  the  output  tape  '.completes  the  process. 

The  coding  structure  of  ICE-III  makes  possible  the  use  of  this 
"slimming  up"  procedure.  This  procedure  permits  the  life-cycle  cost 
estimate  to  include  as  many  levels  as  desired,  since  regardless  of  the 
length  of  the  input  only  nine  levels  of  data  need  to  be  stored  in 
memory  at  any  one  time. 
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6.  DISCUSSION  OF  THREE  PROGRAM  SEQUENCE 


A  sequence  of  three  programs  is  used  to  determine  the  time- 
phased  life-cycle  cost  estimates.  The  three  programs  are: 

•  Data  check  and  Input  Tape  Preparation  Program; 

•  Main  Processing  Program;  and 

•  The  Output  Preparation  and  Summarization  Program. 

These  programs  and  their  functions  are  described  below  in  more  detail. 

6 . 1  Data  Check  and  Input  Tape  Preparation  Program. 

The  functions  of  the  first  program  are: 

•  To  prepare  a  tape  of  data  which  will  be  used  as  input  to 
the  Main  Processing  Program; 

•  To  give  a  print-out  of  data  by  card  type,  which  will 
facilitate  the  location  of  format  errors; 

•  To  check  for  discontinuities  in  the  data.  A  discontinuity 
is  a  difference  of  two  or  more  degree  levels  between  any 
level  and  its  initial  sublevel  (e.g.,  01  and  010101  with  no 
4-digit  code  such  as  0101  in  between) . 

It  is  important  to  note  here  that  the  Msumming  upM  procedure, 
which  is  used  in  the  main  processing  program,  will  work  correctly  only 
if  the  level  data  are  ordered  in  an  increasing  degree  of  complexity. 

The  function  of  ordering  the  level  data,  according  to  the  above  rule,  has 
not  been  included  in  this  program,  due  to  excessive  sorting  time  for  large 
samples  and  the  fact  that  the  data  need  not  be  sorted  for  each  new  run. 
Thus,  an  optional  program  is  provided  which  will  order  the  input  data  and 
provide  a  tape  of  the  ordered  data  for  input  to  the  main  processing 
program. 


Listing  and  flow  charts  of  this  optional  program  are  given  in 
the  appendix.  The  input  for  this  program  assumes  that  the  data  are  in 
the  following  order: 

•  The  three  "L"  cards,  in  order; 

•  The  function  cards,  in  any  order; 

•  A  card  with  10  asterisks  following  the  function  cards; 
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•  The  level  data,  in  any  order;  finally 

•  Two  blank  cards  and  a  PROB  card.* 

6.2  Main  Processing  Program. 

This  program  is  the  heart  of  the  cost-estimating  model;  its 
functions  are : 

•  To  read  the  tape  prepared  by  Input  Tape  Preparation  Program, 
or  the  optional  program; 

•  To  process  these  data  and  obtain  the  life-cycle  cost 
estimates  of  the  weapon/support  system  by  level  and  time 
interval ; 

•  To  write  an  output  tape  containing  the  results  of  the  above 
levels,  which  will  be  used  as  input  to  the  Report 
Preparation  and  Summarization  Program. 

6 . 3  Report  Preparation  and  Summarization  Program. 

This  is  the  final  program  in  the  sequence;  its  functions  are: 

•  To  read  as  input,  the  tape  prepared  by  the  Main  Processing 
Program; 

•  To  sort  this  tape,  by  level,  into  order  of  increasing 
complexity  of  level  number; 

•  To  give  standard  output  of  the  calculated  costs  by  level 
number; 

•  To  provide  an  option  for  obtaining  three  standard  summaries 
of  the  cost  data  by: 

••  major  cost  category; 

••  appropriation;  and 

••  cost  category  by  appropriation. 

It  should  be  noted  here  that  the  three  program  sequence  will 
function  correctly  as  long  as  the  input  cards  are  in  correct  form  even 
though  specific  level  numbers  are  not  associated  with  ICE-III.  However, 


* 

If  the  level  data  are  not  correctly  ordered  the  optional  sorting 
program  must  be  used. 
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if  this  is  the  case,  the  three  standard  summaries,  described  above,  will 
not  preform  correctly,  since  in  the  ICE  format  cost  categories  and 
appropriations  have  been  given  specific  codes. 

Figure  3,  gives  a  flow  diagram  of  the  three  program  sequence, 
and  the  necessary  steps  in  obtaining  a  life  cycle  cost  estimate. 


7.  DESCRIPTION  OF  INPUT  CARDS 

There  are  nine  card  types  which  may  be  used  as  input  to  the 
life-cycle  cost-estimating  model.  Card  types  A,  B,  C,  D,  E,  G,  and  H 
are  used  to  associate  data  with  a  specific  level  of  data  or  group  of 
levels.  The  L  cards  carry  information  pertinent  to  the  introduction  of 
the  entire  study.  The  F  cards  enter  function  data,  which  can  be  used 
by  any  data  level  or  group  of  levels. 

7.1  Summarization  of  Card  Types  and  Their  Uses. 

•  "A"  used  to  associate  a  name  with  each  specific  level 
number.  A  card  type  "A"  must  be  present  with  each  distinct 
level  number.  Two  cards  can  be  used.* * 

•  "B"  used  to  enter  cost  or  quantity  data  into  the  data 
tables.  Each  set  can  contain  as  many  as  five  cards  depend¬ 
ing  upon  which  of  three  input  forms  is  used.** 

•  "C"  used  to  reference  "functions"  which  are  to  be  associated 
with  the  specific  level  number.  From  one  to  five  function 
tables  can  be  referenced  on  each  card,  and  as  many  as  nine 
"C"  cards  can  be  used  at  each  level. 

•  "D"  used  to  enter  data  into  a  special  function  table  FN98. 
The  data  are  stored  here  until  a  new  set  of  "D"  cards  are 
encountered.  Each  set  contains  from  one  to  five  cards 
depending  upon  which  of  the  three  formats  are  used.  (Later 
references  on  "C"  card  to  FN98  will  recall  the  values 
previously  defined  in  "D"  cards). 


★ 

The  second  "A"  card  would  be  a  continuation  of  the  first  "A"  card. 

★  ★ 

There  are  three  distinct  forms  of  "B"  and  "D"  cards.  The  forms  may 
be  combined  when  defining  data  for  a  specific  data  level. 
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Figure  3  Pictoral  Description  of  3  Program  Sequence. 
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•  "E"  used  to  form  a  new  special  function,  FN99,  through 

the  addition  of  one  through  five  separate  function  tables. 
These  "E"  cards  should  be  used  only  when  the  derived 
function  table  is  related  to  only  one  specific  level  number. 
Data  are  retained  in  this  table  until  another  set  of  "E" 
cards  is  encountered.  (Later  references  to  FN99  on  "C" 
card  will  recall  values  defined  in  "E"  cards) . 

e  "F"  used  to  generate  from  1  to  97  different  function 

tables.  Each  set  contains  from  one  to  five  cards  depending 
upon  which  format  is  used.  "F"  cards  are  not  read  or 
associated  with  any  specific  level  number,  and  can  be 
referenced  through  "C"  cards  from  any  level. 

e  "GIJH"  must  be  used  together.  "G"  cards  contain  cost  data 
and  "H"  cards  contain  quantity  data.  "G"  cards  are  used  to 
generate  cost  data  based  upon  learning  curve  calculations. 

e  "L"  used  as  report  header  and  beginning  and < ending  year  of 
the  study.  There  are  three  distinct  "L"  cards  used  with 
each  run,  and  are  located  at  the  very  beginning  of  the  data 
deck. 

The  general  format  of  data  cards  A,  B,  C,  D,  E,  and  G  is: 

Card  Variable 


Columns 

Variable 

Format 

Description 

1-18 

Level  1(1) 
1-1,8 

1811 

The  level  number  associated 
with  this  particular  set  of 
data 

19 

CARDNO 

11 

The  number  of  the  card  which 
is  associated  with  this  type 
for  this  level  number 

20 

Blank 

IX 

21 

Card 

A1 

The  letter  designating  the 
card  type 

22-80 

---- 

Variable 

These  columns  contain  the  data, 
in  different  formats  for 

distinct  card  types,  which 
are  to  be  used  in  the  calcu¬ 
lation  of  this  level  cost 
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7.2  Description  of  Input  Cards  by  Type. 


As  was  stated  previously  columns  1  through  21  contain  the  same 
data  for  card  types  A,  B,  C,  D,  E,  G,  and  H.  These  cards  vary  only  in 
Columns  22-80,  and  so  the  following  description  of  the  card  types  will  be 
concerned  only  with  these  columns. 

"A"  Data  Card  Format:  Used  to  associate  a  name  or  description 
with  each  specific  level  of  data.  At  least  one  "A"  card  must  be  used 
with  each  level,  while  up  to  two  are  allowed. 

Card  Description 

Columns  Variable  Format  of  Variable 

22-80  INAMEl(I),  A9,5A10  Contains  the  name  of  this 

1=1,6  data  level 

If  two  "A"  cards  are  used  the  only  difference  would  appear  in 
Column  19,  where  a  2  would  be  placed  to  represent  that  this  was  the  second 
"A"  card  used  for  this  level,  and  Columns  22-80  would  contain  INAMEl(I), 
1=7,12. 


"B"  Data  Card  Format:  Used  to  enter  data  into  tables  which 
will  be  used  in  the  calculation  of  the  cost  pertinent  to  this  level  of 
data.  The  three  types  of  "B"  cards  discussed  below,  enter  the  data 
into  an  array,  where  each  member  is  associated  with  a  given  interval 
during  the  time  frame  of  the  study. 

B-TYPE  I:  If  the  data  for  this  level  are  constant  for  each 
time  increment  over  the  period  of  the  study,  the  B-TYPE  I  card  should 


be  used. 

Card 

Columns 

Variable 

Format 

Description 
of  Variable 

22-24 

Constant 

A3 

The  word  "ALL" 

25 

Blank 

IX 

26-32 

Value 

F7.3 

The  value  here  will  be 
placed  in  the  BDATA  array 
for  all  time  increments 
during  the  study 

B-TYPE  II:  If  the  data  to  be  entered  through  the  "B"  card 
are  constant  over  some  period  during  time  frame  of  the  study,  then  a 
B-TYPE  II  card  would  be  used  to  enter  the  data  for  this  period. 
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Card 

Description 

Columns 

Variable 

Format 

of  Variable 

22-25 

Constant 

A4 

The  word  "FROM" 

26 

Blank 

IX 

27-28 

IYR1 

12 

Initial  year  of  period 
during  which  the  data  are 

constant 

29 

Blank 

IX 

30-31 

Constant 

A2 

The  word  "TO" 

32 

Blank 

IX 

33-34 

IYR2 

12 

Final  year  of  period  during 
which  the  data  are  constant 

35 

Blank 

IX 

36-42 

Value 

F7.3 

Value  which  is  to  be 
entered  into  the  data 
tables  from  IYR1  to  IYR2 

B-TYPE  III:  If  the 

data  to  be 

entered  with  "B"  cards  vary 

from  year 

to  year,  then  TYPE 

III  would  be 

used. 

Card 

Description 

Columns 

Variable 

Format 

of  Variable 

22-23 

IYR(l) 

12 

Period  during  study  which 
will  have  the  value  given  in 
the  next  field 

24-30 

Value (1) 

F7.3 

Value  to  be  placed  in  B 
table  for  IYR(l) 

31-32 

IYR(2) 

12 

Same  definitions  apply  for 

33-39 

Value (2) 

F7.3 

IYR(I) ,  Value (I) ,  1=1,6 

40-41 

IYR(3) 

12 

42-48 

Value (3) 

F7.3 

49-50 

IYR(4) 

12 

51-57 

Value (4) 

F7.3 
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58-59 

IYR(5) 

12 

60-66 

Value (5) 

F7.3 

67-68 

IYR(6) 

12 

69-75 

Value (6) 

F7.3 

76-80 

Blank 

12 

"C"  Data  Card  Format:  Used  to  reference  functions  containing 
factors  to  be  used  in  the  calculation  of  cost  data  for  this  level.  From 
one  to  five  functions  may  be  referenced  on  each  "C"  card  and  a  maximum  of 
nine  "C"  cards  may  be  used  with  any  one  level. 


Card 

Columns 

Variable 

Format 

Description 
of  Variable 

22-23 

Constant 

A2 

The  letters  "FN" 

24-25 

IFN(l) 

12 

Number  of  the  1st  function 
which  is  referenced 

26-27 

Constant 

A2 

The  letters  "FN" 

28-29 

IFN  (2) 

12 

Number  of  the  2nd  function 
which  is  referenced 

30-31 

Constant 

A2 

The  letters  "FN" 

32-33 

IFN (3) 

12 

Number  of  the  3rd  function 
which  is  referenced 

34-35 

Constant 

A2 

The  letters  "FN" 

36-37 

IFN  (4) 

12 

Number  of  the  4th  function 
which  is  referenced 

38-39 

Constant 

A2 

The  letters  "FN" 

40-41 

IFN  (5) 

12 

Number  of  the  5th  function 
which  is  referenced 

"D"  Data  Cards:  Used  to  enter  either  cost  or  quantity  data, 
used  in  the  calculation  of  costs  pertaining  to  a  specific  level.  There 
are  three  types  of  "D"  cards,  and  they  have  exactly  the  same  format  and 
interpretation  as  the  three  types  of  "B"  cards.*  Therefore,  they  will  not 
be  given  here.  "D"  cards  also  have  another  particular  use.  The  data 

* 

With  the 

obvious  exception 

that  a  "D"  would 

appear  in  Column  21. 
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entered  by  these  cards  are  placed  in  special  function  98,  and  will  be 
stored  there  until  another  set  of  "D"  cards  is  encountered.  Hence, 
these  data  can  be  used  many  times  subsequent  to  their  appearance  on  "D" 
cards  by  referencing  function  98  on  either  "C"  or  "E"  cards. 

"E"  Data  Cards:  Have  a  use  which  is  similar  to  that  of  "C" 
cards.  Recall  that  "C"  cards  reference  functions  which  are  used  as 
factors  in  the  calculation  of  a  particular  level's  cost  data.  "E"  cards 
also  reference  functions.  The  functions  referenced  are  added  by  year  or 
time  increment  and  in  essence  form  a  new  function.  This  new  function  is 
then  used  as  a  factor  in  the  calculation  of  costs  for  the  level. 

Card  Description 


Columns 

Variable 

Format 

of  Variable 

22-23 

Constant 

A2 

The  letters  MFNM 

24-25 

IFN(l) 

12 

Number  of  1st  function 
referenced 

26-27 

Constant 

A2 

The  letters  MFN" 

28-29 

IFN (2) 

12 

Number  of  2nd  function 
referenced 

30-31 

Constant 

A2 

The  letters  "FN" 

32-33 

IFN  (3) 

12 

Number  of  3rd  function 
referenced 

34-35 

Constant 

A2 

The  letters  MFNM 

36-37 

IFN (4) 

12 

Number  of  4th  function 
referenced 

38-39 

Constant 

A2 

The  letters  MFNM 

40-41 

IFN (5) 

12 

Number  of  5th  function 
referenced 

"E"  data  calculated  in  this  manner  are  stored  in  special 
function  99  until  a  new  set  of  "E"  cards  are  encountered.  Thus,  the 
data  can  be  entered  once  then  recalled  by  referencing  function  99  on  a 
"C"  card. 


"F"  Data  Cards  Format:  Used  to  enter  data  into  the  function 
Tables  1  through  97.  (Recall  that  "D"  and  "E"  cards  enter  data  into 
functions  98  and  99  respectively.  Data  input  on  "F"  cards  should  not 
be  placed  in  either  special  function  98  or  99.)  The  data  entered  into 
the  function  tables  are  not  associated  with  a  specific  level  but  can  be 
used  as  factors  in  calculating  the  costs  pertaining  to  any  level. 
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Three  types  of  "F"  cards  are  used  to  enter  data  into  function 
Tables  1  through  97.  The  uses  of  the  three  types  of  cards  are  the  same 
as  those  for  "B"  and  "D"  cards. 

TYPE  1  is  used  when  all  values  in  the  function  are  to  be  a 
single  constant. 


Card 

Columns 

Variable 

Format 

Description 
of  Variable 

1-14 

Blank 

14X 

15-16 

Constant 

A2 

The  letters  "FN" 

17-18 

Function 

Number 

12 

Number  of  function  into 
which  we  are  placing  data 

19 

CARDNO 

11 

Number  of  card  pertaining 
to  this  function 

20 

Blank 

IX 

21 

Constant 

A1 

The  letter  "F" 

22-24 

Constant 

A3 

The  word  "ALL" 

25 

Blank 

IX 

26-32 

Value 

F7.3 

Value  to  be  placed  in  all 
entries  of  this  function 

over  some 

Columns  1- 

TYPE  2  is  used  when  the  value  in 
interval  of  years,  less  than  the 
-21  are  the  same  as  Type  1. 

the  function  remains  constant 
entire  period  of  the  study. 

Card 

Columns 

Variable 

Format 

Description 
of  Variable 

22-25 

Constant 

A4 

The  word  "FROM" 

26 

Blank 

IX 

27-28 

IYR1 

12 

Beginning  FY  of  period 

29 

Blank 

IX 

30-31 

Constant 

A2 

The  word  "TO" 

32 

Blank 

IX 

33-34 

IYR2 

12 

Final  FY  of  period 

35 

Blank 

IX 

36-42 

Value 

F7.3 

Value  to  be  placed  in  the 
function  table  for  this 
period 

43-80 

Blank 

38X 

24 


TYPE  3  is  used  when  the  values  to  be  placed  in  the  function 
table  vary  from  year  to  year.  Columns  1-21  are  the  same  as  Type  1. 

Card 

Columns 

Variable 

Format 

Description 
of  Variable 

22-23 

IYR(l) 

12 

Period  during  study  for 
which  this  function's 
value  will  be  specified  in 
next  field 

24-30 

VAL(l) 

F7.3 

Value  given  this  function 
for  period  given  in  pre¬ 
ceding  field 

31-32 

33-39 

IYR(2) 

VAL (2) 

12 

F7.3 

Definitions  same  as  above 

40-41 

IYR(3) 

12 

42-48 

VAL (3) 

F7.3 

49-50 

IYR(4) 

12 

51-57 

VAL (4) 

F7.3 

58-59 

IYR(5) 

12 

60-66 

VAL (5) 

F7.3 

67-68 

IYR(6) 

12 

69-75 

VAL (6) 

F7.3 

"G§H"  Data  Card  Format:  Used  to  generate  cost  and  quantity 
data  for  a  specific  level.  An  "H"  card  must  be  used  each  time  a  "G" 
card  is  used.  The  "G"  card  contains  the  cost  data,  and  references  a 
function  which  contains  the  quantity  to  be  procured  or  developed.  The 
"H"  card  contains  further  specifications  concerning  quantities.  These 
variables  are  then  used  in  "learning  curve"  calculations  to  determine 
the  costs  for  this  level  during  each  time  increment.  The  uses  of  the 
"G"  and  "H"  cards  will  be  explained  in  more  detail  in  a  later  section. 

Card 

Columns 

Variable 

Format 

Description 
of  Variable 

22-30 

A 

F9.0 

First  unit  cost  for  the 
item  described  in  this 
level 
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31-36 

B 

F6.6 

The  value  of  B  in  learning 
curve  formula 

37-38* 

BPR 

12 

Slope,  which  is  equal  to 
-  lo«  B/-log  2 

39-44 

IQTY1 

16 

Starting  quantity  for  the 
calculations 

45-46 

FN 

A2 

Constant 

47-48 

NUM 

12 

The  number  of  the  function 
which  contains  the  quantity 
data  pertinent  to  this 
level 

"H"  Card 

Format 

Card 

Columns 

Variable 

Format 

Description 
of  Variable 

22-25 

UNIT (1) 

14 

The  number  of  the  unit  which 
describes  the  beginning  of 
the  first  quantity  category 

26-31 

QTY(l) 

F6.3 

The  factor  to  be  multiplied 
by  each  member  of  the  above 
referenced  function,  which 
lies  between  UNIT(l)  and 

UNIT (2) 

32-35 

UNIT (2) 

14 

The  number  of  the  unit  which 
describes  the  end  of  the 
first  quantity  category  and 
the  beginning  of  the  second 
quantity  category 

36-41 

QTY(2) 

F6.3 

The  factor  to  be  multiplied 
by  each  member  of  the  above 
referenced  function,  which 
lies  between  UNIT(2)  and 

UNIT (3) 

42-45 

UNIT(3) 

14 

Same  definition  as  previous 
unit  variables 

*  . 
This  entry  is  superfluous  but  was  included  due  to  its  presence  in 

MICOMS  card  description. 
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46-51 

QTY(3) 

F6.3 

Same  definition  as 
QTY  variables 

previous 

52-55 

UNIT(4) 

14 

Same  definition  as 
UNIT  variables 

previous 

56-61 

QTY (4) 

F6.3 

Same  definition  as 
QTY  variables 

previous 

62-65 

UNIT (5) 

14 

Same  definition  as 
UNIT  variables 

previous 

66-71 

QTY (5) 

F6.3 

Same  definition  as 
QTY  variables 

previous 

"L" 

beginning  and 
and  all  must  1 

Data  Card  Format:  Contains  header  information  and 
ending  fiscal  years  of  the  study.  There  are  three 
be  present  for  each  model  run. 

the 
cards , 

L-l 

: 

Card 

Columns 

Variable 

Format 

Description 
of  Variable 

1-16 

Blank 

16X 

17-18 

Constant 

A2 

The  letters  "LA" 

19 

Constant 

11 

The  number  "1" 

20 

Blank 

IX 

21 

Constant 

A1 

The  letter  "L" 

22-80 

Report  Name 

A59 

L-2 

• 

1-16 

Blank 

16X 

17-18 

Constant 

A2 

The  letters  "LA" 

19 

Constant 

11 

The  number  "2" 

20 

Blank 

IX 

21 

Constant 

A1 

The  letter  "L" 

22-50 

Report  Name 

A29 

Continued 
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51-52 

Day 

12 

Day  of  Run 

53 

Blank 

IX 

54-62 

Month 

A9 

Month  of  Run 

63 

Blank 

IX 

64-67 

Year 

14 

Year  of  Run 

68-80 

Blank 

12X 

L-3: 

1-16 

Blank 

16X 

17-18 

Constant 

A2 

The  letters  "LA" 

19 

Constant 

11 

The  number  "3" 

20 

Blank 

IX 

21 

Constant 

A1 

The  letter  "L" 

22-23 

IYR1 

12 

Beginning  Fiscal  Year 
Study 

24-25 

Blank 

2X 

26-27 

IYR2 

12 

Ending  Fiscal  Year  of 
Study 

7.3 

Further  Explanation 

of  "G" 

and  "H"  Cards. 

The  "G"  and  "H"  cards,  as  stated  previously,  are  used  to 
determine  the  costs  associated  with  a  level  of  data  through  learning 
curve  calculations.  Some  additional  explanation  is  necessary  here  to 
make  clear  the  use  of  these  card  types. 

The  learning  curve  function  has  the  form: 

Y  =  AX"B 

A  =  first  unit  cost  of  the  item  being  produced 

B  *  slope  of  learning  curve 

X  =  quantity  being  produced 

Y  =  cost  of  the  Xth  item 
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Cost/Unit 


INTEGRAL  APPROXIMATION  OF 
LEARNING  CURVE  CALCULATION 


The  total  cost  of  units  through  X0  would  be 


Total  Cost  =  y~l  AX 
X=X, 


-B 


However,  for  large  numbers  of  items  (N>_30) ,  LICEM  uses  the 
integral  approximation: 


Total  Cost 


■? 


AX 


-B 


V1 


ft  M-b; 
i  -  b  r  fXj-i 


The  function  referenced  in  Columns  47  and  48  of  the  "G"  card 
contains  the  procurement  quantities  for  the  time  increments  during  the 
study.  These  quantities  coupled  with  the  quantity  to  start  (Columns 
39-44)  are  sufficient  to  calculate  the  costs  associated  with  each  time 
period. 


If  one  is  calculating  costs  associated  with  the  unit  quantities 
referenced  in  the  function  in  Columns  47  and  48,  the  "H"  card  must  still 
be  used,  but  it  is  really  of  little  value.  For  example,  suppose  that 
function  5  is  referenced,  and  contains  the  number  of  airframes  to  be 
procured  during  the  time  increments  of  the  study,  the  "H"  card  for  this 
example  would  then  be: 
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COL  22-25 


0001 


COL  26-31 


000001 


This  can  be  interpreted  as  the  procurement  of  one  airframe  for  each 
total  missile  to  be  procured. 


Suppose  now  that  in  another  level  the  costs  of  the  fins  to  be 
placed  upon  each  missile  are  to  be  calculated  and  that  there  are  to  be 
four  fins/missile.  In  addition,  five  extra  fins  are  required  for 
testing.  Function  5  would  still  be  referenced  in  this  level,  since 
this  function  contains  the  number  of  missiles  to  be  procured  during  each 
period.  (The  other  entries  on  the  "G"  card  would  be  changed  to  reflect 
the  different  learning  curve  parameters  associated  with  the  production 
of  missile  fins.)  The  "H"  card  in  this  example  would  be: 


COL 

22-25 

0001  i 

COL 

26-31 

000005* 

COL 

32-35 

0006  1 

COL 

36-41 

000004) 

Units  1-5  have  5 
fins/missile 


Units  of  more  than  5  have 
4  fins/missile 


These  data  have  the  following  interpretation: 

Associated  with  units  1-5  of  the  missiles  referenced  in  func¬ 
tion  5,  we  are  using  five  fins/missile.  (These  five  extra  fins  are  those 
desired  for  testing.)  Missile  units  above  5,  have  four  fins/missile 
associated  with  them. 

The  use  of  the  "H"  card  is  now  obvious,  because  through  its 
use  the  quantities  of  each  part  of  the  work  breakdown  structure  need  not 
be  stored  in  the  function  tables.  As  long  as  the  number  of  each 
component  used  in  a  single  unit  is  known,  only  the  number  of  complete 
units  needed  must  be  stored  in  the  function  tables.  This  greatly 
reduces  the  amount  of  data  which  needs  to  be  stored  in  core  during  a 
case . 


Having  described  the  cards  and  their  functions  we  are  now 
ready  to  describe  the  order  of  input  cards  for  this  model.  Figure  4 
pictorially  describes  the  necessary  order  of  input.  Note  the  order: 

1st  -  3  "L"  cards 

2nd  -  All  function  card  data 

3rd  -  Card  with  10*' s  in  1st  10  columns  (used  to  separate 
function  cards  from  level  data) 
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PROB  CARD  SIGNIFIES  THE  END  OF  DATA  FOR 
THIS  RUN  (A  BLANK  CARD  SHOULD  FOLLOW 
THE  LEVEL  DATA  &  SEPARATE  IT  FROM  THE 
PROB  CARD) 


CARDS  FOR  EACH  LEVEL  OF  DATA 
INCLUDED  IN  THE  STUDY  THIS 
DATA  INCLUDES  A,  B,  C,  D,  E,  G, 
AND  H  CARDS. 


SEPARATOR  CARD  SIG¬ 
NIFIES  THE  END  OF 
FUNCTION  DATA  & 
THE  BEGINNING  OF 
LEVEL  DATA. 


Figure  4  Order  of  Input  Cards 
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4th 


-  Level  data,  where  all  levels  are  ordered  in  increasing 
level  complexity. 

5th  -  Blank  card 

6th  -  Prob  Card  (signifies  end  of  problem) 


8.  MISCELLANEOUS  POINTS 

Several  additional  points,  not  logically  falling  within  the 
classification  of  any  of  the  other  sections,  should  be  discussed.  These 
points,  to  be  discussed  in  the  following  sections,  are: 

•  Ease  of  parameter  change  for  sensitivity  analysis 

•  Omission  of  original  data  sort 

•  Special  use  of  "A"  and  "B"  cards 

•  "B-D-F"  card  options 

•  Use  of  previously  made  decks  for  MICOM's  COBOL  version 

•  Summary  of  card  combinations  and  their  uses 

•  Tape  requirements  for  three  program  sequence 

•  Description  of  tapes  produced  by  the  input  tape  preparation 
and  the  main  processing  programs 

•  Compilation  and  running  times  for  the  three  program  sequence 
8.1  Ease  of  Parameter  Change  for  Sensitivity  Analysis. 

The  formats  in  which  the  data  are  entered  into  this  life-cycle 
cost  model  are  especially  suited  to  the  performance  of  sensitivity 
analyses.  If  the  parameter  to  be  investigated  is  stored  in  the  function 
tables,  a  change  of  at  most  five  cards  is  all  that  is  necessary  to 
change  the  variable  value. 

Also,  if  the  variable  of  interest  lies  in  only  one  portion  of 
the  life-cycle  cost  estimate,  say  the  operating  cost  category,  it  is  not 
necessary  to  recompute  the  costs  for  the  other  cost  categories.  This 
can  prove  to  be  a  time  and  money  saving  feature  on  extensive  studies. 
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8.2  Omission  of  Original  Data  Sort. 


In  the  original  COBOL  version  of  this  LCCE  model,  tasks  are 
sorted  into  three  programs.  In  the  initial  program,  similar  to  our 
Input  Tape  Preparation  Program,  a  sorting  routine  ensures  that  the 
data  are  sorted  and  input  in  correct  order. 


In  this  FORTRAN  version  of  the  model,  users  are  provided  with 
two  options.  These  are: 

a.  Data  Check  and  Tape  Input  Preparation  Program  assumes  that 
the  data  are  in  correct  order  with  three  main  functions; 

•  to  prepare  a  tape  containing  input  for  LCCE  model. 

•  to  list  data,  by  card  type,  so  that  a  visual  check  can  be 
made  for  format  errors . 

•  to  check  for  discontinuities  in  the  data. 

b.  Sorter  assumes  that  the  data  deck  is  out  of  order;  this 
program  has  two  main  functions; 

•  to  put  the  data  into  correct  order  for  input  to  LCCE  model. 

•  to  prepare  a  tape  for  input  to  LCCE  model. 


These  options  eliminate  the  need  for  a  complete  sorting  each 
time  a  deck  is  run.  The  sorting  procedure  is  time  consuming  and 
unnecessary  when  a  data  deck  is  to  be  run  more  than  once. 


8.3  Special  Use  of  "A"  and  "B"  Cards. 

It  may  be  that  the  costs  associated  with  some  level  of  data 
may  already  be  known  and  no  calculation  is  necessary  to  determine  them. 
In  this  case  there  are  only  two  card  types  necessary  for  the  cost 
description  of  this  level.  These  are: 


•  "A"  -  Card(s) :  to  describe  the  level. 

•  "B"  -  Card(s) :  to  enter  costs  for  this  level. 


When  this  configuration  of  cards  is  used  to  describe  a  level 
of  data,  the  model  assumes  that  the  costs  entered  on  the  "B"  cards  are 
entered  in  millions,  even  though  the  same  format  (F7.3)  is  used  to  enter 
the  number. 


As  an  example,  if  a  "B"  card  was  to  be  used  to  enter  already 
known  costs,  then  $1,500,000.00  would  be  entered  as  bbbl500,  the  decimal 
point  would  be  assumed  to  lie  between  the  1  and  5. 
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8.4  "B-D-F"  Cards  (Options). 


Note  that  there  are  three  ways  that  data  can  be  entered  through 
the  "B,"  "D,"  and  "F"  cards.  These  ways  can  be  classified  as  "ALL" 
cards,  "INTERVAL"  cards,  and  "YEAR  by  YEAR"  cards.  "ALL"  cards  place  the 
specified  value  into  each  year  of  the  period.  "INTERVAL"  cards  place  the 
specified  value  into  each  year  between  and  including  the  specified  end¬ 
points,  and  the  "YEAR  by  YEAR"  cards  specify  a  new  value  for  each  year 

The  order  in  which  these  cards  appear  is  not  commutative.  That 
is,  the  same  cards  in  different  orders  will  not  give  the  same  results, 
as  is  seen  in  the  examples  below.  (The  example  uses  "F"  cards,  but  the 
same  caution  applies  to  "D"  and  "B"  cards.) 

Example  1* 

FN011  FALL- 1000 
FN012  F732000 

This  combination  will  have  the  following  results: 

70  71  72  73  74  75 

1000  1000  1000  2000  1000  1000 

If  the  two  cards  were  reversed  however,  the  results  would  be: 

70  71  72  73  74  75 

1000  1000  1000  1000  1000  1000 

This  result  caused  by  the  fact  that  the  "ALL"  card  overlays 
its  value  on  top  of  the  2000  entered  by  the  previous  card. 

Consequently,  there  are  many  ways  of  entering  a  data  set  into 
storage  correctly,  but  one  must  keep  in  mind  that  later  cards  overlay 
their  values  on  those  previously  defined. 

8.5  Use  of  Previously  Made  Data  Decks  for  MICOM's  Model. 

Data  decks  which  have  been  prepared  for  use  in  MICOM's  life- 
cycle  cost  model,  although  written  in  COBOL,  can  be  used  as  input  to 
LICEM  with  the  addition  of  one  card  in  the  data  deck.  This  additional 
card  contains  10  asterisks  in  the  first  10  columns.  It  is  inserted  to 
separate  function  cards  from  level  data. 


* 

Assume  time  frame  is  1970-1975. 
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8.6  Summary  of  Card  Combinations  and  Their  Uses. 


In  the  following  discussion,  there  is  always  an  "A"  card 
associated  with  the  level  under  consideration: 


Cards  Present 

Resulting  Computations 

"B"  only 

No  computations  made;  assumed  values  on  "B" 
cards  are  expressed  in  MILLIONS  of  dollars. 

"B$C"  or  "C&D" 

Multiplies  those  values  in  the  "B"  or  "D" 
cards  by  the  values  present  in  the  function (s) 
referenced  on  the  "C"  cards. 

"B5C5D" 

Multiplies  the  values  on  "B"  cards  by  func¬ 
tion  (s)  referenced  on  "C"  card,  and  multiplies 
that  product  by  values  present  on  "D"  card. 

"BSE"  or  "DSE" 

Multiplies  the  values  in  "B"  or  "D"  cards  by 
the  sum  of  the  functions  (>2)  referenced  on 
the  "E"  card. 

"BSDSE" 

Multiplies  the  values  in  the  "B"  and  "D"  cards 
and  multiplies  this  product  by  the  sum  of  the 
functions  referenced  on  the  "E"  cards . 

"GSH" 

Used  to  make  standard  learning  curve  calcula¬ 
tions  of  form 

h 

C  =  X/  AX_B 

x=x2 

"GSHSB"  or  "C" 
or  "D"  or  "E" 

Makes  standard  learning  curve  calculations  and 
multiplies  results  by  data  in  "B"  or  "D"  or  by 
functions  referenced  in  "C"  or  by  the  sum  of 
the  functions  referenced  in  "E." 

This  list  is  not  meant 

to  be  exhaustive  but  it  does  represent  the  basic 

card  combinations  used  for  obtaining  costs  in  a  given  level. 

Two  Binary  Coded  Decimal  tapes  (BCD)  are  needed  for  the  opera¬ 
tion  of  the  three  program  sequence.  One  tape  is  mounted  on  tape  unit  1 
and  the  other  on  tape  unit  2.  Following  is  a  schematic  diagram  of  how 
these  tapes  are  utilized. 
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Tape  Unit  1 


Figure  5  Tape  Requirements  for  Three  Program  Sequence 


8.7  Description  of  Tapes  Produced. 

1 .  The  tape  produced  by  the  Input  Tape  preparation  program 
consists  of  80  column  images  of  the  card  deck  input.  The  arrangement 
and  format  is  as  follows: 

Format 

8A10 
8A10 
8A10 


8A10 


A10,  70X 


8A10 


L  1 

L  2  Cards 
L  3 


Block  of 
Function 
Cards 


10*' s  8  70 
Blanks 


Level  Data 
Including 
All  Card 
Types 


80  Blanks 
Marks  End 
of  Data 
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This  tape  is  produced 
on  unit  1. 


2.  The  tape  produced  by  the  Main  Processing  Program  is  in 
the  following  form: 


Format 


IYR1 ,  IYR2 


(REPNAM(I) ,1=1,9) 
NDAY,  MONTH,  NYR 


ISTORE (J , K) , K=1 , 2 
(LEVEL  NUMBER) 
ISTORE ( J, K) ,K=4, 12 
I DENT.  OF  LEVEL) 


STORE (J , K) , K= 1 , 30 
COSTS  OF  JTH  LEVEL) 


LAST  TWO  RECORDS 
ARE  REPEATED  FOR 
EACH  LEVEL  OF  DATA 


BLANK  RECORD 


END  FILE 


14,  2X,  14 

A9 ,  5A10 ,  A9 ,  2A10 ,  5X 
12,  2X ,  A9 ,  14 

A10,  A8 ,  A9 ,  5A10 ,  A9,  2A10 

(10F12 . 3/10F12 . 3/ 10F12 . 3) 


8.8  Compilation  Time. 

The  compilation  and  running  time  for  the  three  program  sequence 
are  given  below:  (run  times  are  for  example  problem) . 


Compilation  Time 
(min) 

Run  Time 
(min) 

Storage  Requirements 

Sorter 

.47 

.20 

12K 

LCCE 

.90 

1.03 

10K 

Report 

Generator 

.47 

.21 

12K 

Having  all  the  concepts  necessary  to  use  the  model,  we  are  now 
ready  to  work  through  an  example  from  input  preparation  to  a  description 
of  the  final  output. 
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9 .  EXAMPLE  PROBLEM 


In  this  section,  the  time-phased,  life-cycle  cost  estimates  for 
the  fictional  X-4210  vehicle  will  be  derived. 

The  FORTRAN  coding  sheets  used  to  obtain  the  necessary  input 
cards  to  LICEM  are  shown  in  Figures  6,  7,  and  8.  Note  again  the  order 
of  cards,  "L"  cards  followed  directly  by  the  "F"  (function)  cards,  a 
card  with  10  asterisks,  data  by  level  (sorted  into  order  of  increasing 
complexity),  a  blank  card,  and  finally  a  PROB  card,  which  signifies 
the  end  of  data  deck. 

First,  note  in  our  data  set  that  each  level  begins  its  data 
with  an  "A"  card,  and  that  some  levels  have  only  "A"  cards.  When  this 
is  the  case,  the  particular  level,  having  only  "A"  card  identification, 
will  be  the  sum  of  the  costs  associated  with  all  its  sub-levels. 

The  normal  type  of  calculation  that  will  be  required  by  the 
data  will  involve  "B,"  "C,"  and  "D"  cards.  When  combinations  of  "B" 
and  "C"  or  "C"  and  "D"  cards  are  used,  the  data  (cost  or  quantity)  pre¬ 
sent  on  the  "B"  (or  "D")  card  is  multiplied  by  the  data  stored  in  the 
functions  referenced  by  the  "C"  cards.  Note  that  it  is  permissible  to 
reference  a  function  more  than  once  in  a  "C"  card.  This  has  the  effect 
of  taking  that  function  to  a  power.  (See  levels  010103010101  and 
010103010102  for  examples  of  these  calculations) . 

Another  commonly  used  calculation  involves  the  "B"  or  "D"  and 
"E"  cards.  In  this  case  the  data  (cost  or  quantity)  stored  in  the  "B" 
or  "D"  card  are  multiplied  by  the  sum  of  functions  (by  year)  referenced 
in  the  "E"  card.  As  in  the  "C"  card,  a  function  can  be  referenced  more 
than  once  on  an  "E"  card.  This  has  the  effect  of  multiplying  the  values 
in  the  function  by  the  number  of  times  it  is  referenced.  (See  level 
010202  for  example.) 

The  last  standard  calculation  involves  the  "G"  and  "H"  cards. 
These  are  used  to  perform  learning  curve  calculations.  (See  level  010305 
for  example.)  It  is  understood  that  the  function  referenced  in  the 
"G"  card  contains  the  quantities  for  the  learning  curve  calculations. 

This  example  can  also  help  to  make  clear  the  definition  of  a 
discontinuity  in  the  data.  As  defined  earlier,  a  discontinuity  in  the 
data  is  a  difference  of  two  or  more  level  degrees  between  a  level  and 
one  of  its  initial  sublevels.  To  paraphrase  the  definition,  a  level 
which  contains  N  pairs  of  digits  cannot  be  followed  immediately  by  a 
sublevel  containing  N  +  2  pairs  of  digits.  This  input  error  will  cause 
the  "summing  up"  procedure  in  the  main  processing  program  to  operate 
incorrectly,  and  erroneous  results  will  be  given. 


38 


If  discontinuities  do  exist  in  the  data,  SORTER  will  print  out 
a  statement  indicating  the  level  numbers  between  which  the  discontinuity 
exists . 


9 . 1  Preparation  of  Input. 

Notice  in  this  example  that  the  most  detailed  level  of  data 
contains  only  12  digits.  This  level  of  detail  was  sufficient  to 
illustrate  the  important  functions  of  the  program.  The  factors  govern¬ 
ing  the  level  of  detail  for  a  user  are  amount  of  detail  available  and 
amount  desired.  Note  that  the  last  three  pairs  of  digits  are  reserved 
for  the  work  breakdown  structure  of  the  item  being  costed,  and  are  not 
used  in  this  example. 

9 . 2  Output  -  Data  Check  and  Input  Tape  Preparation  Program  (SORTER) . 

Having  described  the  input  for  our  example  case  we  can  now 
discuss  the  output  as  obtained  from  the  three  program  sequence. 

The  objective  of  the  print  out  from  this  sorter  program  is  to 
help  the  user  spot  errors  in  his  input  data  which  could  cause  errors  in 
the  main  processing  model  execution.  Except  for  the  "L"  cards,  all 
other  card  types  are  listed  on  separate  pages  at  the  top  of  which  is  an 
image  of  a  correctly  formated  card,  to  help  spot  errors.  Notice  in 
cards  "B"  and  "D"  and  "F"  where  numbers  are  to  be  entered  three  decimal 
places  are  assumed  and  no  decimal  point  need  be  placed  in  the  field. 
Figures  9  through  23  contain  the  output  from  program  SORTER. 

9 . 3  Output  From  Report  Generator  Program. 

Note  that  the  only  output  of  the  Main  Processing  Program  is 
the  tape  it  produces  for  the  Report  Generator  and  Summarization 
Program. 


The  output  of  the  Report  Generator  and  Summarization  Program 
is  divided  into  two  portions: 

a.  (Listing  of  levels  input  and  their  associated  costs  over 
each  period  during  the  time  frame  of  the  study.)  The  levels  are  in  the 
order  in  which  they  were  input  to  the  model.  (Increasing  order  of 
Complexity) .  In  addition  to  the  cost  by  year,  the  total  cost  attributable 
to  the  level  is  printed  on  the  right  hand  side  of  the  output. 

b.  Optional  Summaries  (obtained  by  placing  a  1  in  column  1  of 
the  first  data  card  for  the  Report  Generator  and  Summarization  Program) . 

(1)  Summary  by  Major  Cost  Category,  where  these  cost 
categories  include  Development,  Investment  Recurring  and  Non  Recurring, 
and  Operating  (Figures  24-25) . 
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(2)  Summary  by  Appropriation,  where  possible  appropriations 
are  RDT§E,  PEMA,  0§MA,  MPA,  MCA,  ASF,  FHMA.  For  this  summary  all  levels, 
which  are  included  in  any  appropriation,  are  summed,  (regardless  of  the 
cost  category  in  which  they  fall)  to  obtain  the  total  for  that  appropria¬ 
tion.  If  any  appropriation  has  a  total  of  0.0,  no  print  out  is  given  for 
that  appropriation. 

From  our  example  the  RDT$E  present  in  the  summary  by  appropria¬ 
tion  is  merely  level  010103010101,  since  it  is  the  only  RDT$E  level 
present . 


The  PEMA  entry  in  this  summary  however  is  the  sum  of  the  costs 
of  levels  010103010102  and  01040101010102.  At  the  top  of  the  page, 
under  SUMMARY  by  APPROPRIATION,  is  given  the  total  of  the  RDT§E  and  the 
PEMA  costs  are  below  (Figure  26) . 

(3)  Finally  the  summary  of  all  appropriations  within  each 
cost  category  is  given  (Figure  27) . 

It  is  possible  that  users  might  desire  other  summaries.  For 
this  reason  the  format  of  the  output  tape  produced  by  the  main  processing 
program  is  given  in  the  Miscellaneous  Points  Section.  Once  this  format 
is  known  the  extraction  of  any  summary  desired  is  relatively  easy. 
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Figure  6  Sample  Problem  Input 
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Figure  7  Sample  Problem  Input 
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Figure  8  Sample  Problem  Input 
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Figure  9  L-Cards  From  Sample 
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SAMPLE  CASE-TIME  PHASED  LIFE  CYCF  COST  ESTIMATING  MOOEL 

COST  DATA  IN  MILLIONS 
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INVESTMENT  RECURRING  0103 

1.094  1*136  1.332  1.294  1,274  1.260  1.250 

1.159  1,167  1,316  1.286  1.269  1.256  1.247 

1,107  1,356  1.304  1.280  1.264  1.253  24,904 


SAMPLE  CASF-TIHE  PHASED  LIFE  CYCF  COST  ESTIMATING  MODEL 

cost  data  in  millions 
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Figure  25  Standard  Cost  by  level 
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SUMMARY  BY  APPROPRIATION 
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summary-cost  category  by  appropriation 
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Next  page  is  blank. 


APPENDIX 


FLOW  CHARTS  AND  PROGRAM  LISTINGS 


This  appendix  contains  flow  charts  and  listings  of  each  program 
in  the  sequence  (see  Figures  A.l  and  A. 2). 
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READ  IN 
'DATA -PLACE 
.DIRECTLY  ON, 
TAPE 


I 


i'read  in  good 

COPIES  OF  15 
DISTINCT  DATA 
CARD  TAPES 


PLACE  THESE 
DATA  IN  ICARD 
ARRAY 
ICARD 

(I,  1,  K),K  =  1,  11 


I 


READ 

IN  & 

PRINT  OUT 

THE 

THREE 

\\|^  n 

CARDS 

READ  IN 
FUNCTION 
CARD 


DETERMINE 
WHICH  TYPE 
OF  FUNCTION 
FORM  THIS  IS 
&  PLACE  IN 
CATAGORY 


YES 


BEGIN  TO 
READ  IN  DATA 
BY  LEVEL 


SORT  THESE  DATA 
INTO  THEIR 
APPROPIATE 
ARRAYS  BY  TYPE 
OF  CARD 


r  PRINT  OUT  DATA 
ICARDS  -SEPARATED 
BY  TYPE 


CHECK  FOR 
DISCON  - 
TINUITIES 
OF  DATA 


I 


GIVE  PRINT  THAT 
INDICATES  DIS¬ 
CONTINUITIES  IF 
THEY  EXIST 


GIVE  LISTING 
OF  COMPLETE 
DATA  DECK 


RETURN 

Figure  A*1  Flow  Chart  for  the  Data  Check  and  Imput  Tape 
Preparation  Program 
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Figure  A-2  Chart  I  -  Flow  Charts  for  Optional  Program 
Which  Sorts  Data  into  Correct  Order. 

THIS  PROGRAM  ASSUMES  THE  FOLLOWING  ORDER  OF  THE 
IMPUT  CARDS : 

1.  THE  3  L  CARDS  (IN  ORDER) 

2.  THE  FUNCTION  CARDS  (NOT  NECESSARILY  IN  ORDER) 

3.  A  CARD  WITH  10  ASTERISKS 

4.  THE  LEVEL  DATA  (NOT  NECESSARILY  IN  ORDER) 

5.  A  BLANK  CARD 
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Figure  A*  2  Chart  H 


58 


A  listing  of  the  Data  Check  and  Input  Tape 
(Program  1  in  the  Sequence)  Follows.  (See  Figures  A. 3 


Preparation  Program 
through  A. 7). 
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CB  17, .29 


FEB,  7,72  BRLFSC2  FORTRAN, 

*  SA11WC  392  107  4210  COPES  PROGRAM  1 

dimension  icardi  u,  ioo,fl)  ,  i tempi  8)  ,ncount<  u)  ,  izi  la) ,  i  zi  no) 
DIMENSION  I  TEMPI ( 8 1 , PLANK  !  8) 

1  FORMAT! A  10, A8, A1 ,A2, Al, A8,5A10) 

2  FORMAT ( 8A10  > 

3  format ( 70X, a i o ) 

4  FORMAT (2A10,2A1,5A10,A8) 

5  FORMAT ( 18A1 ) 

<S  FORMAT  (  IX,  Al) 

7  FORMAT ( 1H1 ) 

R  FORMAT  ( 2 A  1 > 

9  FORMAT! 1H-, 'A  DISCONTINUITY  EXISTS  BETwEEN  LE VELS ’ , 3X, 1 8 A 1 , 3X , ' AND 
1 ' ,3X, 1RA1 ) 

DATA! ICARDI 1, 1,K),K-1, 4)  /'  FNXXl  ' , ' F ALL  XXXXX ' , 

l'XX  •/ 

DAT A ( ICARD(2,1,K),K-1,5)  /'  ','  FNXXl  ','FFROM  FY  T', 

1 '0  FY  XXXXX', 'XX  '/ 

DaTA(ICARD(3,1,K),K*2,8)  /'  FNXXl  ' ,  ' FFYXXXXXXX ' , ' FYXXXXXXXX ' , 
1'  YXXXXXXXFY', 'XXXXXXXFYX' , 'XXXXXXFYXX' , 'XXXXX  '/ 

DAtA(ICARD(4,1,K),K»1,8)  /' —LEVEL  N ',' UMBER - 1  *,'A - ', 

1' - LEVEL  '  ,  '  I  DENT  I  F  I  CA  '  ,  '  TI  ON  NAME-',' - ',' . '/ 

DATA!  ICARD15,  1,K),K»1,4)  /'—LEVEL  N', 'UMBER— 1  ' ,  '  BALL  XXXXX', 

l'XX  '/ 

DATA( ICARD(6,1,K),K-1,5)  /• - LEVEL  N', 'UMBER - 1  *,'BFROM  FY  T', 

1  •  o  FY  XXXXX'/'XX  */ 

DATA!  ICARD(7,1,K),K«1,R)  /'—LEVEL  N', 'UMBER— 1  •  ,  '  BF  YXXXXXXX  •  , 

1'FYXXXXXXXF', 'YXXXXXXXFY* , 'XXXXXXXFYX' , 'XXXXXXFYxX', 'XXXXX  '/ 

DAtA(  ICARD(«,  1,K),K»1,5)  /  '—LEVEL  N',  'UMBER—  1  • ,  '  CFnXXFNXXF  • , 
l'NXXFNXXFNX', 'X  '/ 

DATA! ICARD(9,1,K),K«1,4)  /'---LEVEL  N', 'UMBER-— 1  • , ' DALL  XXXXX', 

l'XX  '/ 

DATA(  ICARD!  10,  1 ,  K  ) ,  K- 1 , 5 )  /  • - LEVEL  N',  'UMBER—  1  *,'DFR0M  FY  T', 

l'O  FY  XXXXX' 9  1 XX  '/ 

DATA!  ICARD!  It,  1,K),K»1,8)/'— LEVEL  N'  ,  'UMBER - 1  '  ,' DF  YXXXXXXX  '  , 

1 'FYxxXXXXXF', 'YXXXXXXXFY' , 'XXXXXXXFYX' , 'XXXXXXFYxX' , 'XXXXX  '/ 

DAtA( ICARD( 12, 1,K),K»1,5)/' - LEVEL  N', 'UMBER - 1  • , ' EFnXXFnXXF • , 

l'NXXFNXXFNX', 'X  '/ 

DATA!  ICARD!  13,  1 ,  K  ) ,  K»1 , 5  >  / ' - LEVEL  N', 'UMBER— i  ' ,  '  OlUNI  TCOST  '  , 

l'BVALUEXXOT', 'YSTRFNXX  '/ 

DATA!  ICARD!  14,  1 ,  K  )  ,  K»1 , 8  )  /  '  —  LE  VEL  N'  ,  '  UMBER— 1  •  ,  '  HlUN-lXXXX  '  , 
l'XXl)N-2XXXX', 'XXUN-3XXXX' , 'XXUN-4XXXX' , 'XXUN-5XXXX' , 'XX  '/ 


DATA  BLANK/10H  / 

DATA  LANK/10H  / 

DATA  ITAR/10H**********/ 
DATA  STAR/10H**********/ 
REWIND  1 
DO  12  1-1,8 


PLANK! I ) -PLANK 
12  CONTINUE 

DO  10  1-1,14 
NCOUNT ( I ) »  1 
10  CONTINUE 
DO  15  1-1,3 

READ  2,  ( ITEMP! J), J-1,8) 
WRITE  (1,2)  ( ITEMP! J), J-1,8) 
15  CONTINUE 

DO  17  1-1,300 

READ  2,  ( ITEMP! J) , J-1,8) 
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o  o 


IF  (ITEMP(1),EQ.ITAR)  GOTO  16 
WRITE  ( 1 r 2 )  ( ITEMPC J) , J-l , 6) 

17  CONTINUE 

16  WRITE  (1,2)  STAR/  ( PLANK ( J > , J-l r 7) 
I8TART-3-M 

DO  20  K I -J  START , 90000 
READ  2,(1  TEMP ( I ) ,  I  ■  1 , 8 ) 

IF  (ITEMP(1),E0.  LANK)  GOTO  25 
WRITE  (1,2)  (ITEMP(I), 1-1/6) 

20  CONTINUE 
C*  *  *  *  <* 

ALL  DATA  NOW  WRITTEN  ON  UNIT  1 


25  WR I TE ( 1 , 2 )  (PLANK(I), 1-1,8) 

END  FILE  1 
REWIND  1 
DO  30  KL-1,3 

READ ( 1 , 2 ) ( ITEMP( I ) , I- 1 ,8 ) 

PRINT  2,  ( ITEHP(I), 1-1,8) 

30  CONTINUE 

DO  100  KL- 1,300 

READ ( 1, 2) ( I TEMP ( I ) , 1-1, 6) 

IF  ( ITEMPC 1 >  #E0*ITAR)  GOTO  110 
DEC0DE(B0,6, I TEMP ( 3 ) )  CHK 
IF  (CHKfEQ,lHA)  GOTO  80 
IF  (CHK.EQ.1HF)  GOTO  50 
NCOUNT ( 3 )— NCOUNT ( 3) *1 
N3-NC0UNT ( 3 ) 

DO  40  1-1,8 
I  CARD ( 3, N3, 1 )-ITEHP( I ) 

40  CONTINUE 
GOTO  100 

50  NCOUNT ( 2 ) -NCOUNT ( 2) +1 
N2-NC0UNT (2) 

DO  60  1-1,6 

I  CARD ( 2, N2, I )  —  ITEMP( I ) 

60  CONTINUE 
GOTO  100 

80  NCOUNTt 1 )-NCOUNT( 1 )  +1 
N 1 -NCOUNT ( 1 ) 

DO  90  1-1,8 
I  CARD ( 1,N1, 1 )-ITEHP( I ) 

90  CONTINUE 
100  CONTINUE 
110  CONTINUE 

READ(1,2)(ITEMP(I), 1-1,8) 

DECODE  (60/5,1  TEMP ( 1 ) )  ( 1 2 ( I ) , I -1 , 1 8 ) 
CALL  DEGREEdZ, LEV, 1/18) 

DO  500  JKL-1, 50000 
DECODE  (80/6, ITEMP(3) )  CHK1/CHK2 
130  IF(CHK1,EQ,1HA)  0OTO  150 
IF(CHKi.EQ,lHB)  GOTO  200 
IF ( CHKl f EQ( 1HC )  GOTO  250 
IF(CHKUEQ.IHD)  GOTO  300 
IF(CHKl,EQ,lHE)  GOTO  350 
IF(CHKi,EQ«lHG)  GOTO  400 
IF (CHKl'EQ* 1HH)  GOTO  450 
150  NCOUNT ( 4 ) -NCOUNT (4)41 
N4-NC0UNT ( 4 ) 

DO  160  1-1,6 
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I  CARD ( 4r N4, 1 ) -ITEMP  < I ) 

160  continue: 

GOTO  480 

200  IF(CHK2tE0#lHA>  GOTO  210 
IF(CHK2,EQ, 1HF)  GO  TO  220 
NC0UNT(7)*NC0UNT(7)*1 
N7*NC01JNT<7) 

DO  205  1*1,8 
ICARD(7,N7, 1 J-ITEMP  C I ) 

205  CONTINUE 
GOTO  480 

210  NCOUNT (5)*NC01JNT(5)M 
N5«NC0UNT ( 5 ) 

DO  215  1*1/8 
I  CARD ( 5,  N5, 1 1 « I  TEMP ( I ) 

215  CONTINUE 
GOTO  480 

220  NC0UNT(6)*NC0UNT(6)+1 
N6*NC01JNT<6) 

DO  225  1*1/8 
I  CARD ( 6, N6, I ) - 1  TEMP  C I  ) 

225  CONTINUE 
GOTO  480 

250  NCOUNT  (8)*NC0IJNT(8)  +  1 
N8*NC0UNT(8) 

DO  295  1*1/8 
I  CARD ( 8, N8, 1 > - ITEMP  C I ) 

295  CONTINUE 
GOTO  480 

300  IF(CHK2,EG,1HA>  GOTO  310 
IF<CHK2,ECI,1HF)  GOTO  320 
NCOUNT ( 1 1 )*NCOUNT ( 1 1 )+l 
N11*NC0UNT<11> 

DO  305  1*1/8 
I  CARD ( 1 1 / N1 1 / I > ■ I  TEMP ( I ) 
305  CONTINUE 
GOTO  480 

310  NCOUNT ( 9 )*NCOUNT 1 9 1 ♦ 1 
N9*NC0IJNT(9) 

DO  315  1*1/8 
ICARD(9,N9/  n-ITEMP(I  ) 

315  CONTINUE 
GOTO  480 

320  NCOUNT ( 1 0 1 "NCOUNT ( 10  J  +1 
N 1 0*NC0UNT ( 10) 

DO  345  1*1/8 
I  CARD ( 1 0/ N 1 0/ 1 > ■ ITEMP  < I ) 
345  CONTINUE 
GOTO  480 

350  NCOUNT ( 12) "NCOUNT ( 12) ♦! 
N12-NC0UNT< 12) 

DO  395  1*1/8 
I  CARD ( 12/N12/ I ) - ITEMP ( I ) 
395  CONTINUE 
GOTO  480 

400  NCOUNT! 13)*NC0UNT( 13)+1 
N13-NC01JNT!  13) 

DO  405  1*1/8 
I  CARD ( 13/N13/ 1 ) - ITEMP ( I ) 
405  CONTINUE 


GOTO  480 

450  NCOUNT ( 14 ) -NCOUNT ( 14 )  ♦  1 
N14-NC0UNT< 14) 

DO  475  IM,8 
ICARD(14,N14,  I)MTEMP(I) 

475  CONTINUE 

480  READ  (1/2)  ( I  TEMPI ( I ) , IM , 8 ) 

IF  ( I  TEMPI ( 1 ) .EQ.LANK )  GOTO  510 
DECODE  (80/5, I  TEMPI ( 1) ) ( I Z 1 ( I)/ 1*1/ 18) 

DO  485  KM, 18 

IF ( IZ ( K ) .NE.IZl ( K > >GOTO  488 

485  CONTINUE 

486  DO  4&7  KM, 8 

I  TEMP ( K ) M  TEMPI ( K ) 

487  CONTINUE 
GOTO  500 

488  CALL  DEGREEUZl, LEVI, 1,18) 

NDIF-LEV1-LEV 
IF(NDIF;GE,2)G0T0  490 

491  DO  489  KM ,  & 

ITEMP(K)MTEMPKK) 

489  CONTINUE 
LEV-LEV1 
GO  TO  500 

490  PRINT  9,((IZ(I),IM,1&)/(IZ1(I)/IM/16)) 

GO  TO  491 

500  CONTINUE 
510  DO  600  JKM,  14 
PRINT  7 

IKNOW-NCOUNT( JK) 

DO  550  JOB-1, IKNOW 

PRINT  2,  (  ICARD(  JK,  J0B,M),MM,8) 

550  CONTINUE 
600  CONTINUE 
REWIND  1 
STOP 
END 

SUBROUTINE  DEGREE  ( I  ARRAY, NDEG, I  ST ART, I  STOP ) 
DIMENSION  I  ARRAY ( 18 ) 

DO  10  I ■ I  START , I  STOP 
IF  ( I  ARRAY ( I ) «EQ« 1H  )  GO  TO  15 
10  CONTINUE 
NDEG-9 
GO  TO  20 
15  NDEG- ( I-l )/2 
20  CONTINUE 
RETURN 
END 

*  L I  ST ( STOP ) 
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100 


Figure  A-3  Chart  I 
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Figure  A-3  Chart  II 
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AT  THIS  POINT  THE  FINDAT 
ARRAY  CONTAINS  THE  COSTS 
FOR  THE  LEVEL  WE  HAVE 

JUST  PROCESSED.  MUST  NOW 
DETERMINE  HOW  THIS  DATA 
SHOULD  BE  STORED  BY 
COMPARISON  WITH  THE  DEGREE 
OF  THE  NEXT  LEVEL  OF  DATA 

ADD  UP  AND  WRITE  ON 
TAPE  ALL  LEVELS  OF 
GREATER  COMPLEXITY 
THAN  LEVI,  WHICH  ARE 
PRESENTLY  BEING  KEPT 
IN  STORE  MATRIX.  THIS 
MUST  BE  DONE  BECAUSE 
DATA  NOW  STORED  IN 
STORE  (LEVI,  K),  K  =  l,  30, 
WILL  BE  REPLACED  BY 
LATEST  LEVEL  OF  PRO¬ 
CESSED  DATA. 


7^ 


'CALL  BACK 
TO  PLACE 

(  1200  ) - ►(  LEVEL  NUMBER 

JNTO  IZ(1)  & 

12(2) 


_ /  \  IN 


FOR  1=1,2 
ISTORE  (LEVI,  I)= 
12(1) 


NEXT 

PAGE 

- ►- 


Figure  A_3  Chart  III 
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Figure  A~3  Chart  1ST 
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BECAUSE  THE  NEXT  LEVEL  IS 
OF  LESS  COMPLEXITY,  WE 
MUST  ADD  UP  &  WRITE  ON 
TAPE  ALL  LEVELS  FROM  LEVI 
TO  LEVI -2  INCLUSIVE 


WE  HAVE  ENTERED 
AND  PROCESSED 
OUR  LAST  LEVEL 
OF  DATA 


I 


SUM  UP  &  WRITE 
ON  TAPE  ALL 
ENTRIES  IN  STORE 
ARRAY  WHICH  HAVE 
NOT  ALREADY  BEEN 
PLACED  ON  TAPE 


FOR  1  =  1,  30 
STORE  (LEV2,  I)  = 
STORE  (LEV2,  I)  + 
STORE  (LEV2-1, 1) 


I 


PRINT  KTAPE 
THE  TOTAL  NUM¬ 
BER  OF  LEVELS 
PLACED  ON 
TAPE 


TRANSFER 
TEMP  2  TO 
TEMPI  ARRAY 


Figure  A-3  Chart  Y 
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A  listing  of  the  Main  Processing  Program  (Program  2  in  the 
sequence)  follows.  (See  Figure  A. 8) 


7,72  BRLFSC2  FORTRAN, 

SA11WC  392  107  4210  COPES  PROGRAM  2 


CB  17, .31 


FEB, 


S 

S 


1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
18 

17 

18 

19 

20 


25 


COMMENT  USES  TAPE  1  AS  INPUT 
COMMENT  USES  TAPE  2  AS  OUTPUT 

DIMENSION  REPNAM(9),FUNnAT<a>,VAL(6),TEMPim,TEMP2(a>*!FN(5>, 
AD AT  A (  12) ,BDATA(3fl) , CDATA(30) , dDAT a ( 30 ) , ED AT  A <  30 ) ,  GD AT A  (  30  ) 
DIMENSION  DATLEVI30) ,12(2) , F 1 NDAT ( 30 >, LEVEL  1 ( 18) 

DIMENSION  LEVEL2I 18) , I YR(5> 

DIMENSION  FUNC ( 99/ 30 ) /  I NAHE1 (  12)  ,  FN<5>,  OT Y {  10),  NUNIT(6) 
DIMENSION  I  STORE (9, 15) /STORE (9, 30) 

DIMENSION  TTF.MPI24) 

DATA  XER/IOHXXXXX  / 

DATA  BLANK/10H  / 

DATA  LANK/10H  / 

FORMAT  ( 8 A  1 0 ) 

FORMAT  (6X/I2/3X/A1) 

FORMAT  I  lX/6( I2/F7.3) ) 

FORMAT  (5X/F7.3) 

FORMAT  I6X, I2/4X/ 12/ 1X/F7.3) 

FORMAT  ( 21X, A8/5A10) 

FORMAT  I21X/ A9/2A10/ 12/ IX, A9/ IX, 14) 

FORMAT  I  2 1  X , 1 2, 2X, 12) 

FORMAT  ( 18R1,A1,1X/2A1) 

FORMAT  I1X/A9/5A10) 

FORMAT  ( 3X, 12, 4 ( 2X, 12)  ) 

FORMAT  ( 1X,F9,0,F6.6/2X/ I6,2X, 12) 

FORMAT  I  IX, 5( 14/F6.3) ) 

FORMAT  I1H1,AB,5A10,A8,2A10,/,50X,I2,2X,I2,2X,I2) 

FORMAT  (I4/2X/I4) 

FORMAT  I A9,5A10, A9,2A10,5X,I2/ 1X/A9/1X, 14) 

FORMAT  (A10/172X) 

FORMAT  I A10, A8, A9/5A10, A9/2A10) 

FORMAT  I  1  OF  12, 3,/, 1  OF  12, 3, 7/1  OF  12. 3) 

FORMAT  115) 

KTAPE-0 

IB»0SIC»O!*IJ)«0SIE«OSIG»O 
REWIND  1 
REWIND  2 
MI NER-99999 

READ  (1,1)  ( TTEUP ( I ) , I  *  1 , 8 ) 

READ  11,1)  (TTENP( I ), I»9, 16) 

READ  (1,1)  (TTEMP( I ) , I-17/24) 

DECODE (80, 6, TTEMP(  1 )  )  (REPNAMt I >, 1-1/6) 

DECODE ( 60, 7 , TTEHP (  9))  ( ( REPN AM ( I > , I *7, 9 ) , NdA Y , NmOnTH, N YR ) 

DEC0DE(80,8/TTEMP( 17) )  I YRl, IYR2 

WRITE  12/16)  I  I REPNAM( I ) / 1*1/9) / NDAY/NMONTH/NYR) 

DO  25  1*1,9 
I  STORE ( I , 3 ) ■ 1 
CONTINUE 

NTOTYR*! YP2-I YR1+1 


WRITE  (2,15)  IYR1/IYR2 
50  READ  (1,1)  ( FUND AT ( I ) , 1*1 , 8) 

C*****  CARD  WITH  10*  MUST  FOLLOW  FUNCTION  CARDS 

IF  (FUNDAT( 1 ) .EG.10H********** >  SO  TO  90 
DECODE  ( BO , 2, FUNDAT ( 2 ) )  NUM/TYPE 
IF  ( TYPF.,FQ,  1HA  )  GO  TO  70 
IF  ( TYPF.,FQ,  1HF )  GO  TO  80 

DECODE  180/3, FUNDAT ( 3 ) )  <(IYR(I),VAL(I)),I*1,6) 
DO  60  1*1,6 
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IF  (  I  YP  (  I  >  .F-G.lOH  )  GO  TO  50 

J-IYR<  I I-IYR1+1 
FUNC(  NIIM,  .T »  -  V  A  L  (  I  ) 

60  CONTINIiR 
GO  TO  50 

70  DECODE  (60, 4,FUHDAT<3)  )  VALU 
DO  75  I « 1 , NTOT YR 
FUNC <  NUM, I ) ■VALU 
75  CONTINIJF. 

GO  TO  50 

80  DECODE  <  80 , 5, FUNDAT ( 3  )  )  NSTART,  NEND,  VAlU 

1 1- NSTART-IYRl+l 

12- NEND-IYR1+1 
DO  85  1-11,12 
FUNC ( NUM, I > -VALU 

85  CONTINIJF 
GO  TO  50 

C*  *  *  *  * 

C  AT  THIS  POINT  WF  ARE  REApY  TO  bFGIN  PROCESSING  LEVEL  DATA 

C-**** 

90  DO  91  1-1,9 

I  STORE ( I , 3) a  1 

91  CONTINIJF. 

100  READ  (1,11  (TEHn  (  I  )  ,  1-1,8) 

102  decode  ( an, g, tempi ( 1 ) >  t ( level i ( n , i»i, ifl) ,cardno. card,  type > 

CALL  DEGREE  (LEVELl,LFVl, 1,181 
IF  (CARD.rQ, 1HA)  GOTO  105 
IF  ( C ARR.FQ . 1HB )  GOTO  111 
IF  (CARD.FQ.IHC)  GOTO  138 
IF  (CARD.FQ.IHD)  GOTO  160 
IF  ( CARD .FQ. 1  HE )  GOTO  183 
IF  ( CARD.FQ. IMG )  GOTO  200 

C-**** 

C  THIS  SFCTION  PROCESSES  A  TYPE  DATA  CaPDS 

C  CARNO  TFLLS  WHETHER  1ST  OR  2ND  A  CARn  FOR  THIS  LEVEL 

105  I  A -  1 

IF  (CARDNn.EG.lH2)  GO  TO  HO 

DECODE  <80, 10, TEMPI (3)  1  ( I  NAME  1 (  I  >, 1-1,6) 

GO  TO  1000 

110  DECODE  <80,10,  TEMPI (3) )  ( I  NAME  1(11,1-7,12) 

GO  TO  1000 

Q  ★  ★  *  ★  # 

C  THIS  SECTION  PROCESSES  Ft  TYPE  DATA  CaPDS 

C  *  *  *  *  * 

111  1 8-1 

IF  (TYPF.FO. 1HA)  GOTO  120 
IF  ( TYPF.FO. 1HF )  GOTO  130 

DECODE  <80,3, TEMPI ( 3  > )  ( ( I YR < I ) , V AL < I ) ) , I ■ 1 , 6 ) 

DO  112  1-1,6 

IF  ( I YR« I ) .FG.10H  ) GOTO  lOOO 

J  —  I YR ( I ) - 1 YR 1  + 1 
BDATA ( J ) -VAL  < I > 

112  CONTINUE 
GOTO  1000 

120  DECODE  (80, 4, TEMPI (3) )  VALU 
DO  125  I - t , NTOTYR 
BDATAI I ) -VALU 
125  CONTINIJF. 

GOTO  1000 
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130  DECODE  <80, 5, TEMPI (3) )  NST ART, NEND , VALu 
I 1-NSTART-IYR1+1 
I  2-NEND- 1 YR 1  + 1 
DO  135  1*11,12 
BD AT  A ( I )-VALU 
135  CONTINUF 
GOTO  1000 

C  THIS  SECTION  PROCESSES  C  TYPE  DATA  CARDS 

138  IC*1 

DECODE  ( 80 ,11, TEMP  1(3)1  <  I  FN < I) , I *1 , 5 ) 

IF  (CAPDN0.NE.1H1 )  GOTO  141 
I ■ I FN (  1  ) 

DO  140  J-t, NTOTYR 
CDATA( J)-FUNC( I, J) 

140  CONTINUE 
KL-2 

GOTO  142 

141  KL*  1 

142  DO  145  I *KL , 5 

IF  ( I FN ( I ) .EQ.LANK )  GOTO  1000 
J-IFNC I) 

DO  143  K*  t , NTOTYR 
CDATA(K)-CDATA(K)*FUNC( J,K) 

143  CONTINUE 
14$  CONTINUE 
150  GOTO  1000 


C  THIS  SECTION  PROCESSES  D  DATA  CARDS 


160  ID-1 

IF  (CARDN0.NE.1H1  >  GOTO  163 
DO  161  1*1, NTOTYR 
DDATA ( I >-0. 

161  CONTINUE 

163  IF  (TYPF.EQ.1HA)  GOTO  170 
IF  (TYPF.EQ.1HF)  GOTO  180 

DECODE  <80,3, TEMPI ( 3 ) )  I (I YR (I) r VAL< I) ) , I -1 , 6 ) 
DO  162  1-1,6 

IF  ( I YR ( I ) .E0.2H  )  GOTO  1000 

J-I YR( I I-IYR1+1 
DDATAl J)-VAL( I > 

162  CONTINHF 
GOTO  1000 

170  DECODE  (80, 4, TEMPI (3) )  VALU 
DO  172  1-1, NTOTYR 
DD ATA ( I >*VALU 
172  CONTINUE 
GOTO  1000 

180  DECODE  (80, 5, TEMPI (3) )  NST ART, NEND/ VALu 
I  1 -NST  ART- 1 YR 1 +1 
I 2*NEND-I YR 1+1 
DO  182  1*11,12 
DD ATA ( I >-VALU 
182  CONTINUE 
GOTO  1000 

C***** 

C  THIS  SECTION  PROCESSES  E  DATA  CARDS 

C***** 
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183  I E- i 

IF  (CARDNO.NE. 1H1 )  GOTO  184 
HO  186  I  - 1 , NTOTYR 
ED ATA ( I ) - 0  . 

186  CONTINUE 

184  DECODE  (80, 1 1, TEMPI (3) )  ( I FN ( I ) , I ■ 1 , 5 ) 

I  -  I FN ( 1 ) 

DO  185  J- 1 , NTOTYR 
EDATA( J)«FUNC< I, J) 

186  CONTINUE 

DO  1<)0  1-2,5 

IF  < IFN( I ) ,EQ.2H  >  GOTO  1000 
J-  I  FN ( I ) 

DO  188  K-t, NTOTYR 
EDATA(K)-EDATA(K)+FUNC( J,K) 

188  CONTINUE 
190  CONTINUE 
GOTO  1000 


THIS  SECTION  PROCESSES  G  AND  H  DATA  qARDS 


200  I G  *  1 

DECODE  (80, 12,TEMP1(3)  )  A ,  B,  I  <JT  Y  1  ,  I  FN 
READ  (1,1)  (TEMR2( I), 1-1,8) 

DECODE  (80, 13,TFMP2(3> )  ( ( NUNIT ( I ) , QTY < I ) > , I ■ l , 5 ) 
DO  210  1-1,5 

IF  (NUNJT(I).EQ.IOH  >  GOTO  211 

210  CONTINUE 
NALL-5 

NUNIT (6)-10000000000 
GOTO  213 

211  NUNIT ( I >-10000000000 
213  CONTINUE 

NTOT 1 -I QTY 1 
DO  215  1-1,5 
J-  I  ♦  1 

IF  (NT0T1 ,LT.NUHIT( J)  )  GOTO  218 
215  CONTINUE 
218  I  CL- I 

DO  290  IYR-1, NTOTYR 
TRACK-0. 

I  OIJ  AN-FUNC  (  IFN,  I  YR) 

NTOT2-NTOT14-IOUAN 

IF  (NT0T2, EQ, NTOTl >  GOTO  290 

IF  ( NT0T2.GT.NUM I T ( I CL+1 ) )  GOTO  230 

IF  ( IQUAN.LE.30)  GOTO  225 

CALL  WEDDLE ( A, B, NTOTl , NT0T2, QTY ( I  CL > /TOTAL) 
TRACK-TRACK+TOTAL 
GDAT  A ( I YR ) -TRACK 
NTOTl -NT0T2 
GOTO  280 
225  NBOOB-NTOTl 

CALL  SUMUrt A,B,MB00B,NT0T2,QTY( I  CL), TOTAL) 

TRACK -TRACK+TOTAL 
GDAT  A ( I YR ) -TRACK 
GOTO  280 

230  DO  250  JJ-I CL, 5 

IF  (NT0T2,LE.NUHIT( JJ+1 ) )  GOTO  260 
NUMB-NUN IT( JJ+1) -NTOTl 
IF  (NUMB. OT. 30)  GOTO  235 
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NB00B-NTDT1 

CALL  SIJMIJP  <A,B,nBOOB,nUnIT(JJ+1),QTY<  JJ>,  TOTAL) 

TRACK-TRACK+TOTAL 

GOTO  240 

23b  CALL  WEDDLE  ( A, R, NTOT 1 , NUN  I T ( JJ+I) , QT Y < JJ ) , TOT AL > 
TRACK-TRACK+TOTAL 
240  NTOTl-NTOTl+NUMFUl 
250  CONTINUE 
260  NUMB-NT0T2-NT0T1 

IF  (NUMB. GT. 30)  GOTO  265 
NBO0B-NTOT1 

CALL  SIJMUP  (A,B,NB00B,NT0T2,QTY(JJ), TOTAL) 

ICL-JJ 
GOTO  270 

26b  CALL  wFDDI.E  (  A,  B,  NTOT  1 ,  NTOT2,  QTY  (  JJ  ),  TOTAL  ) 

I  CL* J  J 

270  TRACK-TRACK+TOTAL 
GDATA ( I YR ) *TR ACK 
280  NTOT 1*NT0T2+1 
7.00  CONTINUE 
300  CONTINUE 
GOTO  1000 


AT  THIS  POINT  HE  HAVE  PROCESSED  A  CARD, 

IF  THF.  NEXT  DATA  CARD  IS  OF  THE  SAME  LEVEL  WE  WILL  PROCESS  IT, 
IE  NOT  WE  CAN  DETERMINE  ITS  DEGREEf  PERFORM  THE  CALCULATIONS 
ON  THE  PREVIOUS  LEvELf  AND  STOPE  THE  RESULT  ACCORDING  TO  HOH 
THE  LEVEL  DEGREES  COMPARE 

•  •  •  *  * 

1000  RE  AD ( 1 , 1 )  ( TEMP2 ( I ) , I *1 , 8 ) 

DECODE  ( 80 , Q, TEHP2 ( 1 ) )  ( LE VEL2 ( I ) , I ■ 1 f 1 8 ) 

DO  1005  I ■  1 , 1 8 

IF  (LEV ELI ( I ) ,NE.LEVEL2< I ) )  GOTO  1 0 1 0 
100b  CONTINUE 

DO  1006  1*1,8 
TEMPI ( I ) "TEMP2 ( I ) 

100ft  CONTINUE 
GOTO  102 


IF  HERE  WE  GOTO  102  WE  ARE  STILL  PROCESSING  DATA  FROM  THIS  LEVEL 
IF  WE  GO  TO  1010  WE  H A VF  COMPLETED  PROCESSING  CARDS  FROM 
THIS  LFVF.I.  AND  MUST  CALCULATE  THE  COSTS  OF  THIS  LEvEL 

*  *  *  *  * 

1010  call  DEGPFE  (LEVEL2,LEv2, 1, 18) 

IF  (IB.FO.O.AND.IC.EQ.O.AND.ID.EO.O. 

1  AND. IF.EO.O.AUD. IG.EQ.O)  GO  TO  1031 
IF ( IA.EQ.l . AND* IB.EQ. I.ANr, IC.EQ.O.AND, ID.EQ.O.AND.  IE.Eq.O, 

1  AND, IG.EO.O )  GOTO  1034 
DO  1011  I  *  1 , NTOT  YR 
F I NDAT ( I ) ■ 1 , 

1011  CONTINUE 

IF  (IB.F.0,0)  GOTO  1015 

DO  1012  I  *  1 , NTOT YR 

FI NDAT ( I ) -FI NDAT ( I ) *  BOAT  A ( I ) 

1012  CONTINUE 

1015  IF  (IC,EO,0)  GOTO  1020 
DO  1017  I  *  1 , NTOT YR 
FI NDAT ( I >*FINDAT( I >*CDATA(  I  ) 

1017  CONTINUE 

1020  IF  ( I DfEQ,0 )  GOTO  1025 


74 


o  o  o  o  o  o 


DO  1022  I  a  1 , N1 OTYR 
FUNC(98, I ) aDDATA ( I ) 

F I NDAT U ) -FI NDAT (I  )  *  DDAT  A  (  I  ) 
1022  CONTINUE 

1025  IF  (IEfF.O,0)  ROTO  1030 
DO  1027  I-l,NTOTYR 
FUNC(99,I ) »EDAT  A ( I ) 

F I NDAT ( I ) »F I ND AT ( I )*EDATA( I ) 
1027  CONTINUF 

1030  IF  (IGtFQ,0)  GOTO  1035 
DO  1032  I ■ 1 , NTOT YR 
FINDATd  )»FINDAT(  I  )*GDATA(  I  ) 

1032  CONTINUF 
GOTO  1035 

1031  DO  1033  I  a i , NTOT  YR 
F I NDAT ( I ) ■ 0  * 

1033  CONTINUF 
GOTO  1035 

1034  DO  1037  I ■ 1 , NTOTYR 

FINDAT ( I )«BDATA( I )*1000000» 
1037  CONTINUF 

1035  DO  1036  I  a  1 , NTOT YR 

FI NDAT ( I)aFINDAT( I  1/1000000* 

1036  CONTINUF 
I  A  ■  0 

I  B  ■  0 
ICaO 
ID-0 
IE-0 
IG-0 


AT  THIS  POINT  FINDAT  ARRAY  CONTAINS  THE 
DATA  FOP  THF  LEVEL  we  have  .tust  PROCESSED) 

MUST  NOW  determine  HOW  THIS  LEVEL'S  data 
SHOULD  RE  STORED  BY  COMPARISON  WITH  THE 
DEGREE  OF  THE  NEXT  LEVEL, 

DO  1038  1*1, NTOTYR 
BDATA ( I )»0, 

CDATA ( I )»0, 

DD AT  A ( I ) ■ 0 , 

ED AT  A ( I >«0, 

GDATAI I )»0, 

1038  CONTINUF 

IF  (  LE  V2  ,F.Q  ,  0  )  KDONE-1 

IF  ( IST0RF(LEV1,3> ,EQ,1)  GO  TO  1200 

KTAPE-KTAPEtl 

K0P»9-LF.V1  +  1 

DO  1130  I  ■  1  ,  KOP 

J-9- I ♦ 1 

IF  (  ISTORFI J, 3) ,EQ,1)  GO  TO  1130 
DO  1126  Kal,NTOTYR 

STORE! J-1,K)»ST0RE( J-1,K > ♦STORE! J,K) 

1126  CONTINUE 

WRITE  (2,18)  (  (  ISTORE! J,K) ,K-1,2> , ( ISToRE! J,K) ,K»4, 12) > 

WRITE  (2,19)  (STORE! J,K> ,K"1, NTOTYR) 

1  STORE ( J, 3 ) a  1 
1130  CONTINUE 

1200  CALL  PACK(LEVEL1( 1), IZ, 18) 

I  STORE ( LEV  1 , 1 ) - 1 Z (  1 ) 
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I  STORE ( LEV  1 , 2 )  ■  I Z  (  2 ) 

ISTORE(LFV1,3)-0 
DO  1205  1*1,12 
IJK- I +3 

1  STORE ( LEVI, I JKI-INAHE1C 1 ) 

INAME1 ( I ) *L  ANK 
120b  CONTINUE 

DO  120ft  K- 1 , NTOTYR 
STORF (LFV1 ,K)*FINDAT(K) 

120ft  CONTINUF 

IF ( KDONF ,FQ  .  1  )  GOTO  1320 
IF(LEV1,GF.LEV2)  GOTO  1220 
DO  1210  1*1,8 
TEMPI ( I ) -TEMP2 ( I ) 

1210  CONTINIJF 
LEV 1*LEV2 
GOTO  102 

1220  IF(LEV1,GT.LEV2)  GOTO  1250 
KTAPE-KTAPE+1 

WRITE  (2,18)  ( ( I  STORE <  LEVI , K ) , K- 1 , 2) , ( I  STORE ( LEV  1 , K ) , K-4, 1 2 ) ) 
WRITF  (2,19)  (STOPE(LFVl,K), K-l, NTOTYR) 

I  STORE ( LEV  1 , 3) *1 

HINUS*LFVl-l 

DO  1230  1*1, NTOTYR 

STORE  (MINI'S,  1  ) -STORE  (MINUS,  I  )  ♦STORE  (  LEV  1 ,  I  > 

1230  CONTINUF 

DO  1240  1*1,8 
TEMPI  ( I ) -TEMP2 ( I ) 

1240  CONTINIJF 
LEV  1 *LFV2 
GOTO  102 

1250  NDIF-LF.V1-L.EV2 
J-LEV1 

DO  1270  KKK*1,NPIF 
K  T  APE*KT  APE* 1 

WRITE  (2,18)  ( ( ISTORE( J,K) ,K-1,2) , ( IST0RE( J,K ) ,K-4,  12) ) 

WRITF  (2,19)  (ST0RE(J,K),K*1, NTOTYR) 

I  STORE ( J, 3 ) *  1 
DO  1260  K*l, NTOTYR 

STORE ( J-1,K ) -STORE ( J- l, K) ♦STORE ( J,K) 

1260  CONTINIJF 
J-J-l 

1270  CONTINUE 

IF(LEV2.E0,1)  GOTO  1300 
DO  1280  K-l, NTOTYR 

STORE  (LF.V?-1,K) -STORE  (LEV2-1,K)*ST0RE(LEV2,K) 

1280  CONTINUF 
1300  KTAPE-KTAPE+1 

WRITE  (2,18)  (  (  IST0RE(LFV2,K) ,K*1,2) , ( IST0RE( LEV2,K )  ,K-4,  12) ) 

WRITE  (2,19)  (ST0RE(LEV2,K) , K-l, NTOTYR) 

IST0RE(LEV2,3)-1 
DO  1310  1-1,8 
TEMPI ( I ) -TEMP2 ( I ) 

1310  CONTINIJF 
LEV  1 -LEV2 
GOTO  102 
1320  J-LEV1 

DO  1340  I-1,LEV1 
KT  APE-KT  APE  +  1 

WRITE  (2,18)  ( ( IST0RE( J,K ) ,K-l,2) , ( ISToRE( J,  K ) ,  K-4, 12) ) 


76 


WRITE  (2,19)  <STORE< J,K) ,K"i,NTOTYR) 

IF  (J-l,F-OfO>  GO  TO  1340 
DO  1330  K-1,NT0TYR 

STORE< J-1,K > -STORE (J-1,K) ♦STORE <  J,K) 

1330  CONTINUE 

I  STORE ( J, 3 ) ■ 1 
J-J-l 

1340  CONTINUE 

WRITE  (2,17)  XER 
END  FILF  2 
STOP 
END 

SUBROUTINE  WEDDLE  ( A , R, I  1 , 1 2, F ACTOR, TOT AL ) 

BOT-1 .-R 

CON-A/BOT 

F  I  RST" ( FLO AT ( 12) )**( 1*-R) 

SECOND" ( FLOAT ( 1 1 ) ) **( l.-B) 

F0RCE"F I RRT -SECOND 
TOT AL-FORCE* CON* FACTOR 
RETURN 
END 

SUBROUTINE  SIJHIIP  (  A ,  B,  I  1 ,  I  2,  F  ACTOR,  TOT  AL  ) 
TOTAL-O. 

IFUl.EQ.O)  I  1  - 1 
DO  10  1-11,12 
X  I  - 1 

TOTAL"TOTAL+< A*xi**(-R)  ) 

10  CONTINUE 

TOTAL"TOTAL*FACTOR 

RETURN 

END 

SUBROUTINE  DEGREE  ( I  ARRAY, NDEG, I  START, I  STOP ) 
DIMENSION  I  ARRAY ( 25 ) 

DO  10  I ■ I  START , 1  STOP 
IF  ( I  ARRAY ( I ) ,E0. 1H  )  GO  TO  15 
10  CONTINUE 
NDEG"9 
GO  TO  20 
15  NDEG«( I  - 1 ) /2 
20  CONTINUE 
RETURN 
END 

-  L I  ST ( STOP ) 
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Figure  A- 4  Flow  Chart-  Report  Preparation  and  Summarization 
Program 
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A  listing  of  the  Report  Preparation  and  Summarization  Program 
(Program  3  in  the  sequence)  follows. 
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7,7 2  BRLFSC2  FORTRAN, 

SA11WC  392  107  4210  COPES  PROGRAM  3 


CB  17, .32 


FEB, 


S 


1 

2 

A 

4 

5 
ft 
7 

10 
1  1 
14 
lb 
16 

17 

18 

19 

20 
21 
22 

23 

24 
2b 
2ft 

27 

28 

29 

30 

31 

32 

33 

34 

35 
3ft 
37 
3ft 


comment  tape  unit  2  used  for  input 

DIMENSION  m<9) , I  DATE (3) ,nC0UNT<3) ,NYR(30) 
n I  MENS  I  ON  VALUE!  10  0,43), IT (25) , I  NT ( 25) 

DIMENSION  TEMPO (4, 7, 30), TOTAL (4, 7), COScAT (30), APPRO (30) 
DIMENSION  EMT(6) ,GMT( 4) ,HMT(4) ,FMT1 (3) ,GMT1 ( 3) ,HmTK  5) 
FORMAT ( I4,2X, 14) 

FORMAT <5A 10, l OX, 12, IX, A9, IX, 14) 

FORMAT  ( A10, A8,A9,5A10,A9,2A10) 

FORMAT  (  1  OF  12,3,/, 1  OF  12,3,/, 1  OF  12. 3) 

FORMAT ( 1H l, 120X, 'PAGE', IX, 13) 

FORMAT  ( A9/5A10, A9,2A10,  5X, 1 2/ 1 X , a9, 1 X , I  4 ) 

FORMAT  ( 40X, 1  COST  DATA  IN  MILLIONS*) 

FORMAT ( lH0,A9,5A10,55XfAl0fA8) 

FORMAT  ( 1H0, A9,5A10,/, A9, 2 A 1 0, 85X , A  1 0 , A» ) 

FORMAT* 1M0,120X, 'PAGE', IX, 13) 

FORMAT  (ftri0,3) 

FORMAT  ( 2  X , 12, 16) 

FORMAT  ( 1M0,2X, 'DEVELOPMENT' ) 

FORMAT  ( 1U0,2X, ' INVESTMENT-NON  RECURRING') 

FORMAT  (  11(0, 2X,  '  INVESTMENT  RECURRING') 

FORMAT  ( 1H0,2X, 'OPERATING* ) 

FORMAT  ( 1U0, ' SUMMARY  BY  APPROPRIATION') 

FORMAT  (12,16) 

FORMAT  ( 1M0,2X, ' RDT-E' ) 

FORMAT  ( 1H0,2X, ' PEMA ' ) 

FORMAT  ( 1H0,2X, ' 0-MA' > 

FORMAT  ( 1M0,2X, ' MPA' ) 

FORMAT  ( 1H0,2X, ' MCA' ) 

FORMAT  ( 1H0,2X, ' ASF' ) 

FORMAT  ( 1H0,2X, 'FHMA' ) 

FORMAT  ( 1H0, ' SUMMARY  BY  MAJOR  COST  CATEGORY') 

FORMAT  ( 2X , I  2 ) 

FORMAT  ( 1H0, ' SUMMARY-COST  CATEGORY  BY  APPROPRIATION') 


FORMAT  (12) 

FORMAT  (  1 n  A  1  ) 

FORMAT  ( 8 A  1 ) 

FORMAT  (110) 

FORMAT  ( 2X, I  2 ) 

FORMAT  (  I  2, I  ft ) 

integer  FHT,FMT1 ,gmt, gmti,hmt,hmti 

DATA  (FMT ( I ) , I»l,6)  /8H(1H  ,2X,,1H  , 1  OH < ' F Y ' , 1 2/ 6/ 3HX ) , , 1H  , 

1  10HX, 'TOTAL' )/ 

DATA  (GMT < I ) , I«l>4)  /8H(1H  ,3X,/lH  , 10H( ’Ey' , 12, 6f3HX) )/ 
nATA  ( HMT (  I  ) , I  *  1 , 4 )  /8H(1H  ,4X,,1H  / 1 0  H ( ' FY ' , I2/6/3HX)  )/ 

DATA  (FMT1 ( I ) , 1-1,3)  /RHI1H  ,2X,,1h  , 1  OH ( F9, 3, 2X ) ) / 

DATA  ( GMT  1(1), 1*1,3)  /BH ( 1 H  /3X,,1H  , 1  OH ( F9, 3, 2X ) ) / 

DATA  (HMTl U ), 1-1,5)  /ftH<lH  ,4X,,lH  , 1 qH ( F9 . 3 , 2X ) , , 1 H  ,8HX/Fl0.3)/ 


REWIND  2 

READ  (2,6)  ( ( I D( I ) ,1*1,9) , ( I  DATE ( I ), 1*1,3) ) 
READ  (2,1)  I YR1 , I YREND 
NTOTYR-I YREND-IYR1+1 
NM0D»M0D(NT0TYR,3) 

IF  (  NMOD,F.O,0  )GOTO  110 
IF(NMOD,EO,l )GOTO  120 
N  C  0  U  N  T ( 1 ) « ( NTOTYR/3 )  +  1 
NCOUNT (2)*NCOUNT ( 1 ) 

NCOUNT  (  3  )  -  NCOIJNT  (  1  >  - 1 
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ROTO  140 

110  NOD- ( NTOTYR/3 ) 

DO  112  1-1,3 
NCOIJNT  (  I  )  -NOD 

112  continue 

GOTO  140 

120  N COUNT ( 1 )-(NT0TYR/3) +1 
NCOUNT(2)-NCOIJNT(  1 )  - 1 
N COUNT (3)-NC0UNT (2) 

140  LOBO-NCOUfU  ( 1 ) 

MOBO-NCOIINT(2) 

N0B0«NC0UNT(3) 

IF(LOB0,ro,10)  GOTO  141 
FMT(2)-L0G0 

FHT1 ( 2)»L0B0 
GOTO  142 

141  FMT ( 2 ) -  1  12 
FMT 1 ( 2 ) *  1  12 

142  IF (MOBO.EO,  10  )  GOTO  143 
GMT (2I-HOB0 

GMTl ( 2 ) -MOBO 
GOTO  144 

1 43  GMT (21-112 
GMTl (21-1  12 

144  IF(NOBO«EOf 10)  GOTO  146 
HMT (21-N0B0 

HMT 1 ( 2 1 -NO BO 
GOTO  147 

146  HMT (21-1  12 
HMT 1 ( 2 ) -1 1  2 

147  DO  145  I-1,L0B0 

N YR  (  I ) ■ I YP 1 ♦ ( 3* ( I - 1 )  ) 

145  CONTINUE 
L2-L0B0+ 1 

L3-L2  +  NC0IINT(  21-1 
J-  1 

DO  148  1-1.2, L3 

NYR( I )-( I YR1+1 )+( J-l ) *3 

J-J+l 

14ft  CONTINUE 
J-l 

L4-L3+1 

L5-L4+N0B0-1 

DO  150  I  -1.4 ,  L5 

N YR ( I )-( IYR1+2)*! J-l 1*3 

J«J«-1 

150  CONTINUE 

DECODE! 10,36, FMT ( 2 ) )  NFMT 
NSPACF-1 17-2-NFMT*  10 
ENCODE! 10,1, FHT(5) )  NSPACF 
nECODE( 10,36,HMT1  (2) )  MFMT 
HSPACF-l  17-4-MFHTM  1 
ENCODE! 10, 1,HHT1( 4) )  MSPACF 
NFINAL-1  1 +NTOTYR 
DO  200  1-1,1000 

READ  (2,3)  (VALUE! I, J), J-l, 11) 

IF ( VALUE ( 1,1) ,E0.5HXXXXX)G0T0  210 
READ  (2,4)  ( VALUE! I, J), J-12,NFINAL) 
200  CONTINUE 
210  NUM-I-1 
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DO  212  J«l,NUM 

DECODE ( AO, 34, VALUE < J, l))  <IT(K),K-1,10) 

DECODE (BO, 35,VALUE( J,2) )  ( IT ( K ) , K-l 1 , 1 8 ) 

DO  211  K-i,ia 
I  NT ( K  )  - 1 T  ( K  ) 

IF  (  IT(K)#NF..1H0)  GOTO  211 
I  NT ( K ) “0 

211  CONTINUF 

ENCODE (fiO, 34, VALUE! J,42) )  ( I  NT ( K ) , K- 1 , 1 0 ) 
ENCODE<ao,35,VALUE< J , 4  3 ) )  ( I  NT ( K ) , K- 1 1 , 18 ) 

212  CONTINUF 

CALL  ARANGI (VALUE, NUM, 43, 100,43) 

CALL  ARANGI ( VALUE, NUM, 42, 100,43) 

*  *  * 

num  contains  total  number  of  levels  used 

VAI.IIF  NOW  CONTAINS  THE  ORDERED  LEVELS  AND  DATA 

*  *  * 

Inum*i 

DO  300  I  PAGE* 1 , 500 
PRINT  5, I  PAGE 

PRINT  6, ( ( I D ( I ) , 1-1,9) , ( I DATE< I ) , 1-1,3) ) 

PRINT  7 

PRINT  FMT  ,  ( NYP  < I ) , I«1 , LOBO ) 

PRINT  GMT  ,  (NYP( I ) , I*L2,L3) 

PRINT  HUT  ,  (NYR< I ), I«L4,L5) 

I  NIJM9  -  I  NlJM  +  9 
DO  250  I-INUM, IUUMQ 
DO  215  J-Q, 1 1 

I F  (  VALUE.  ( 1 ,  J)  ,NE  •  1H  )  GOTO  225 
215  CONTINUF 

PRINT  10, ( ( VALUE ( I , J), J«3,8) , ( VALUE ( I , j) , J-1,2) ) 

GOTO  230 

225  PRINT  11,  ( ( VALUE( I , J) ,J«3, 11 ) , ( VALUE( I , J) , J-1,2) ) 

230  TO-O. 

DO  231  K0- 1 2, NF I N AL 
TO«TO+VALUE( I,KO) 

231  CONTINUE 

PRINT  FMT1, ( VALUE ( I  ,  K  )  , K-12, NFI NAL, 3) 

PRINT  GMT  1 , ( VALUE ( I , K ) , K« 1 3, NF I N AL,  3 ) 

PRINT  HMT 1 , ( (VALUE(IfK),K-l4,NFINAL/3),T0) 

IF ( I . EQ,NUM ) GOTO  310 
250  CONTINUE 

PRINT  14,1  PAGE 
INIJU-INUHOM 
300  CONTINUE 
310  CONTINUF. 

READ  16,  UOY 
IF  (NOYtF.Q#0>  GOTO  410 


C  CALCULATION  of  summary  1 

C  R Y  COST  CATEGORY 

c  *  *  *  *  * 

I  PAGE- I  PAGE* 1 
PRINT  5,  T PAGE 

PRINT  6, ( ( ID( 1), 1-1,9), ( I  DATE ( I ), 1-1,3) ) 
PRINT  FMT  ,  ( N YR ( I ) , I « 1 , LOBO ) 

PRINT  GMT  ,  (NYR( I), I-L2,L3) 

PRINT  HMT  ,  (NYR( I ), I-L4,L5» 

J-l 

1-1 
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312  IF  (J.F.Q.5)  GOTO  325 

DECODE  <80, 16, VALUE! I, l> )  I  TEST  1 , I TEST2 
IF  !  ITF.ST1.NE.J.0R.ITEST2.NE.0)  OOTO  320 
IER-1 

DO  315  K-12,NFINAL 

TEMPO ( J, 1, IER)-TEMPO( J, 1, I ER ) +V ALUE < I,K) 

I  ER- 1 ER+1 
315  CONTINUE 
J-  J  +  l 
1-1*1 
GOTO  312 
320  1-1*1 

GOTO  312 

c  *  *  *  *  * 

C  TEMPO  ARRAY  CONTAINS  DATA  FOR  1ST  SUMMARY 

C  WILL  CALCULATE  TOTAL  ARRAY  AND  TOTAL  BY 

C  YEAR,  THEN  PRINT  THE  RESULTS. 

C***** 

325  DO  340  1-1,4 

DO  335  K»1,NT0TYR 

TOTAL (If  1 ) -TOTAL ( 1 1  1 ) ♦TEMPO! I ,  If K) 
COSCAT(K)-COSCAT(K)+TEMPO( I,  1,K> 

335  CONTINUE 

340  CONTINUF 
PRINT  30, 

DO  341  I - 1 , NTOT YR 
COSTOT-COSTOT+COSCAT! I ) 

341  CONTINUE 

PRINT  FMT 1 ,  (C0SCAT(K),K-1,NT0TYR,3) 

PRINT  0MT1, (C0SCAT(K),K-2,NT0TYR,3) 

PRINT  HMTl, (C0SCATCK),K-3,NT0TYR,3),C0sT0T 
DO  350  J- l , 4 

GOTO  ( 342f 343, 344, 345) , J 

342  PRINT  17, 

GOTO  346 

343  PRINT  1R, 

GOTO  346 

344  PRINT  19, 

GOTO  346 

345  PRINT  20, 

346  PRINT  FMTl, (TEMPO! J, 1 , K> » K- 1 , NTOT YR f 3) 

PRINT  GMT l, (TEMPO! J, 1 , K > , K-2 , NTOT YR , 3 ) 

PRINT  HMT 1 , ( TEMPO ( J  , 1 , K  » , K-3 , NT OT YR , 3 ) , TOT A L ( J , 1 ) 
350  CONTINUE 


CALCULATION  OF  SUMMARY  2 
BY  APPROPRI ATION 


DO  348  1-1,4 
DO  348  J-1,7 
TOT  AL  (  I ,  J  )  -O', 

DO  348  K-1,30 
TEMPO ( I , J, K ) -0 , 

COSCAT ( K ) -0 , 

34ft  CONTINUF. 

COSTOT-O. 

DO  370  L- 1 , NUM 

DECODE  (80,22,VALUE(L,2) )  I  TEST  1 , I TEST2 
DO  355  J-1,7 

IF  ( lTESTl,EQ.J.ANDtlTEST2tEQ,0)  GOTO  360 
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35b  CONTINUE 
GOTO  370 

360  IER-1 

DO  361  K»12,NFINAL 

APPPOUFR)  -APPRO  <  IER)+VALUE(L,K> 

TEMPO! 1,J, I ER) -TEMPO! I , J, ! ER > +VALUE ( L, K ) 

IER-IER+1 

361  CONTINUE 
370  CONTINUE 

DO  390  J-1,7 
DO  390  K-l, NTOTYR 

TOTAL!  If J) -TOTAL! 1, J>+TEMpO< 1, Jf K) 

390  CONTINUE 

DO  391  K-l, NTOTYR 
APTOT-APTOT+APPRO! K ) 

391  CONTINUE 

IF  !  APTOT.FQ.O. )  GOTO  411 

IPAGE-IPAGE+I 

PRINT  5,  IPAGE 

PRINT  6, ! ! ID! I )f I-lf  9) , ! IDATE! I ) f I-lf 3)  ) 

PRINT  FMT  ,  (NYR< I > , 1-1, LOBO) 

PRINT  GMT  ,  (NYr  «  I  >  f  I-L2,  L3) 

PRINT  HMT  ,  (NYP! I ) , I-L4,L5) 

PRINT  21, 

print  FMT If ( APPTO! K) f K-l , NTOTYRf  3 ) 

PRINT  GMT l, ( APPRO! K) , K-2,  NTOTYR, 3) 

PRINT  HMT 1 , { APPPO! K) , K-3, NTOTYRf  3 ), APTOT 
DO  410  J-1,7 

IF (TOTAL! 1,J) . EQ  .  0 . )  GOTO  410 

GOTO  (402,403, 404, 405, 406, 407, 408), J 

402  PRINT  23, 

GOTO  409 

403  PRINT  24, 

GOTO  409 

404  PRINT  25, 

GOTO  409 

40b  PRINT  26, 

GOTO  409 

406  PRINT  27, 

GOTO  409 

407  PRINT  28, 

GOTO  409 

408  PRINT  29, 

409  PRINT  FMT1, (TEMPO!  1, J,K» ,K-1,MT0TYR,3) 

PRINT  GMT  1 , (TEMPO! 1, J,K> , K-2, NTOTYR, 3 ) 

PRINT  HMT 1 , ( TFMPO ( 1 , J, K ), K-3, NTOTYR, 3 ) , TOTAL! 1, J) 

410  CONTINUE 
*  *  *  *  * 

CALCULATION  OF  SUMMARY  3 

COST  CATEGORY  RY  APPROPRIATION 

★  *  *  ★  * 

411  I  PAGE- 1 P AGE+ 1 
PRINT  5,  IPAGE 

PRINT  6,  ( ( ID! I ) , I-lf 9) , ( I  DATE! I ) f I-lf 3) ) 

DO  420  1-1,4 
DO  420  J-1,7 
TOTAL! I, J)-0. 

DO  420  K-l, NTOTYR 
TEMPO! I, J,K)-0, 

COSCAT ( K ) -0, 
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420  CONTINUE 
COSTOT-O, 

PRINT  FMT,  (  N  YR  (  I) ,  1-1,1.080) 

PRINT  GMT  ,  (NYR(I) ,  I-I.2,L3) 

PRINT  HMT,  (NYR{  I  ) , I-L4,L5> 

HO  430  1-1,4 
00  425  J-1,7 
no  424  l.»l,NUM 

DECODE  (80,37, VALUE <  L, l  M  ITEST1 
DECODE  (80,38, VALUE(L,2) )  1TEST2, ITEST3 

I F  < ITESTi,NE.I,OR,ITEST2*NE.J,OR.ITEST3.NE.O)  GOTO  424 
IER-1 

DO  412  K-t2,NFINAL 

TEMPO ( I, J, IER)-TEMPO( I , J, I ER ) ♦ VALUE < Lr K ) 

IER-IER+1 
412  CONTINUE 

424  CONTINUE 

DO  421  K- t , NTOTYR 

TOTAL ( I, J) -TOTAL ( I , J > *TEMpO < I , J , K ) 
COSCAT(K)-COSCAT(K)+TEMPO( I, J,K) 

421  CONTINUE 

425  CONTINUE 
430  CONTINUE 

DO  422  K-l, NTOTYR 
COSTOT-COSTOT+COSCAT (K ) 

422  CONTINUE 

IF  (COSTOT.EQ.O. )  STOP 
PRINT  32, 

PRINT  FMT 1 ,  (COSCAT(K), K-l, NTOTYR, 3) 

PRINT  GMT  1 ,  (COSCAT(K) ,K-2, NTOTYR, 3) 

PRINT  HMT  t ,  (COSCAT(K) ,K-3, NTOTYR, 3) ,COSTOT 
DO  500  1-1,4 
I  PR  I  NT-0 
DO  480  J-1,7 

IF  (TOTAL! I, J> ,EO.O,  )  GOTO  480 
I F ( IPRINT.EQ.1 )  GOTO  460 
GOTO  ( 456, 457, 458, 459) , I 

456  PRINT  17, 

I  PR  I  NT  -1 
GOTO  460 

457  PRINT  18, 

IPRINT  -1 
GOTO  460 

458  PRINT  19, 

IPRINT  -1 
GOTO  460 

459  PRINT  20, 

IPRINT  -1 

460  GOTO  (461,462,463,464,465,466,467), J 

461  PRINT  23, 

GOTO  470 

462  PRINT  24, 

GOTO  470 

463  PRINT  25, 

GOTO  470 

464  PRINT  26, 

GOTO  470 

465  PRINT  27, 

GOTO  470 

466  PRINT  28, 
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GOTO  470 
467  PRINT  29, 

470  PRINT  FMT I ,  ( TEMPO ( I , J, K ) , K- i , NTOTYR,  3 ) 

PRINT  GMT  1 r  (TEMPO! I, J,K),«»2,NT0TYR,3) 

PRINT  HMTl,  (TEMPO! I, J,K>,K»3,NT0TYR,3),T0TAL(  I,  J> 

480  CONTINUE 
500  CONTINiJF 
STOP 
END 

SUBROUT  INF  ARANGI ( I  A , N, L, NMAX, LM AX > 

0 1  MENS  I  ON  IA(NMAX,LMAX) 

ORPF.RS  LTH  ELEMENTS  OF  ARRAY  IA  FROM  SMALLEST  TO  LARGEST  VALUE 
THE  OTHER  (LMAX-1)  ELEMENTS  OF  ARRAY  IA  ARE  CARRIED  ALONG 
ONE  TO  ONE  AS  THE  LTH  ELEMENT  IS  SEQUENCED 
I F ( N.LE, 1 ) RETURN 
00  40  12-2, N  S  I1-I2-1 
IF  C  I A  (  I1,I..)#LE,IA(I2,U  >GOTO  4  0 
no  10  LL»  1  /  LMAX 

IT-IA(Il,|.L)  S  IA(I1,LL>»IA(I2,LL>  S  I  a  (  I  2»  LL )  ■  I T 
10  CONTINUE 
J2- I  1 

20  I F ( J2 . LF *  1 ) GOTO  40  S  J1-J2-1 
IF  (  IA(J1,|.)«LE#IA(J2»L) )  SOTO  40 
no  30  LL«  1 , 1.MAX 

IT-IA(J1,|.L>  S  I  A  (  J1,LL>»!A(  J2,LL>  *  IA<J2,LL)»IT 
30  CONTINUE  %  J2-J2-1  S  GOTO  20 
40  CONTINUE  %  RETURN  S  END 
L I  ST ( STOP ) 
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